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SECTION  I 


INTRODUCTION 

This  manual  describes  for  the  engineering  user  the  FORTRAN  IV  soft- 
ware used  in  the  Engineering  Flight  Simulation  Validation  study,  contract 
F33615-77-C-2065.  Section  II  includes  an  overview  of  the  programming 
strategy,  and  descriptions  of  each  subroutine.  Section  III  describes 
the  form  of  the  input  and  output  to  this  set  of  programs  and  presents 
sample  input  and  output  files.  Section  IV  provides  a complete  example 
for  a typical  application  of  this  software. 
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SECTION  II 

SOFTWARE  OVERVIEW  AND  SUBROUTINE  DESCRIPTIONS 

The  organization  of  the  programs  in  this  software  package  is  divided 
into  two  sections  — scenario  generation,  and  drive  logic  stimulation/ 
measures  generation.  The  next  sections  discuss  the  FORTRAN  subroutines 
which  implement  the  scenario  generation.  These  are  followed  by  the 
subroutines  for  drive  logic  s imulat ion/ measure s generation. 

MAINLINE 

The  purpose  of  the  scenario  generation  mainline  routine  is  to  set  up 
an  array  of  variables  to  be  plotted  and/or  used  as  inputs  to  the  simulated 
drive  logic.  After  initialization,  the  program  reads  a list  of  the  possible 
variable  names  for  the  particular  problem.  The  program  allows  for  up  to 
150  variable  names  — 50  states  (Z),  50  derivatives  of  the  states  (ZD)  and 
50  auxiliary  variables  (XAUX).  Thus,  a total  of  150  variable  names  must 
be  read  into  the  NVLA  array,  although  many  of  these  names  will  be  blanks. 
Next,  the  number  of  variables  to  be  plotted  is  read  (NV).  If  NV  is  greater 
than  20,  an  error  message  is  printed  and  execution  is  halted.  The  name 
of  each  variable  to  be  plotted  is  then  read,  along  with  its  corresponding 
plotting  scale  factor,  in  units/inch  (KVLA,  SF).  Finally,  the  total  length 
of  the  time  history  (TMAX,  in  seconds)  and  the  time  scale  factor  (TIMSF, 
in  seconds/inch)  are  read.  The  program  then  checks  each  plotted  variable 
name  (KVLA)  against  the  list  of  all  possible  variable  names  (NVLA)  — if 
a plotted  variable  name  is  not  found  in  the  NVLA  list,  an  error  message 
is  printed  and  execution  halts.  If  all  KVLA  names  are  correct,  the  driving 
routine  for  the  scenario  generation  (SENGEN)  is  called. 

ERRMES 

This  routine  is  called  whenever  a fatal  error  is  encountered  during 
execution  of  any  routine  in  the  software  package.  It  prints  out  a particular 
error  message  based  on  the  value  of  IERR.  After  ERRMES  is  called,  execution 
is  halted,  and  the  user  must  correct  the  error  before  attempting  to  rerun 
the  Job. 
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OUTR 


This  subroutine  is  called  during  the  scenario  generation  phase.  Its 
purpose  is  to  buffer  the  data  and  interact  with  the  plotting  subroutine 
GPDPLT . After  initialization,  GPDPLT  is  called  with  10  = 0,  to  plot 
the  axes.  The  data  is  then  stored  in  the  XOUT  array  for  each  variable 
which  was  chosen  to  be  plotted  in  the  mainline  routine.  If  the  maxi mum 
time  has  been  reached  (T  = TMAX)  or  the  10  counter  indicates  that  the 
XOUT  array  is  full  (l00  data  points  for  each  variable)  the  accumulated 
segment  for  each  variable  is  plotted.  Otherwise,  the  routine  returns  to 
its  calling  routine  so  that  the  next  data  points  can  be  generated  and 
buffered. 

GPDPLT 

This  routine  uses  CALCOMP  plotting  routines  to  plot  the  time  histories 
for  the  desired  variables.  The  data  is  transmitted  from  OUTR  through  the 
DATA  array.  On  initial  call,  with  NPTS  = 0,  the  axes  are  plotted,  five 
to  a page.  Thereafter,  the  data  is  plotted  when  GPDPLT  is  called,  in 
groups  of  100  data  points  per  variable.  The  direction  (left  or  right) 
in  which  the  data  is  plotted  is  dependent  upon  the  position  of  the  pen 
so  as  to  minimize  pen  movement.  When  OUTR  sets  the  switch  IEND  in  GPDPLT' s 
argument  list  to  (l  ),  it  signals  to  GPDPLT  that  the  data  in  the  DATA  array 
is  the  last  to  be  plotted.  When  this  data  has  been  plotted  a label  is 
appended  to  the  time  history,  the  axes  are  labeled  with  variable  name 
and  scale  factor,  and  the  time  scale  is  drawn.  Control  returns  to  OUTR. 

For  additional  information  on  the  CALCOMP  subroutines,  refer  to  the  CALCOMP 
Plotter  Manual. 

8ENGEN 

This  subroutine  is  the  driving  routine  for  the  scenario  generation. 

After  initialization,  the  scenario  number  is  read  (see  Section  III  for 
scenario  number  definitions).  Next,  the  number  of  discrete  data  points 
for  the  particular  scenario  is  read,  and  then  the  arrays  of  maximum  times, 
roll  accelerations,  pitch  accelerations  and  derivative  of  inertial  accelera- 
tion which  define  the  inertial  characteristics  of  the  trajectory  are  input 
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(TMAX,  PD,  QP,  VTPP) . The  data  is  generated  at  0.1  second  intervals. 

Within  each  interval,  DATGEN  is  called  to  generate  the  remaining  scenario 
data,  the  data  is  written  to  a file  for  use  as  a forcing  function  by  the 
drive  logic  simulation  routine,  the  units  of  the  data  are  changed  from 
radians  and  g's  to  degrees  and  ft  for  plotting  purposes,  and  OUTR  is 
called  to  buffer  the  data  for  plotting.  This  process  continues  until 
the  maximum  time  (TMAX)  is  reached.  Control  is  returned  to  the  mainline 
routine. 

DATGEN 

This  routine  implements  the  kinematic  equations  which  generate  the 
variables  for  the  scenario.  The  initial  conditions  are  determined,  based 
on  the  scenario  number.  The  states  and  state  derivatives  (Z0,  ZDO)  are 
updated.  Then,  a subroutine  is  called,  based  upon  the  scenario  number, 
which  selects  the  correct  discrete  value  of  roll  acceleration,  pitch 
acceleration  and  derivative  of  inertial  acceleration  (PD,  QP,  VTPP)  based 
on  the  current  time  T.  (in  this  package,  the  subroutines  called  are  HQPT4 
for  the  HQPT4  scenario  and  READAT  otherwise).  The  kinematic  equations  are 
then  computed  to  generate  state  derivatives  (ZD)  as  well  as  ax  and  az 
( 0T ( 1 ) and  0T(2)].  In  addition,  the  variables  (corresponding  to  actual 
flight)  used  as  forcing  functions  by  the  drive  logic  simulation  are  com- 
puted (PA,  QA,  RA,  UPOD,  VPOP,  WPOP) . The  state  derivatives  are  then  "inte- 
grated," using  a second-order  Adams -Bashf or th  numerical  approximation, 

(see  Appendix)  to  arrive  at  the  states  (Z).  Whenever  an  exact  algebraic 
expression  is  available  for  computing  any  state,  it  is  used  following  the 
numerical  approximation  which  it  then  supersedes.  Control  returns  to 
DATGEN. 

READAT 

This  routine  is  used  for  scenario  generation  whenever  the  inputs  to 
the  kinematic  equations  (PD,  QP,  VTDD)  can  be  expressed  as  discrete  values. 
The  purpose  of  this  routine  is  to  compare  the  current  time,  T,  with  the 
current  TMAX  value  — if  T is  less  than  TMAX,  the  current  EP,  QP  and  VTDP 
values  are  used.  Otherwise,  the  counter  INDEX  is  incremented  and  the  next 
discrete  value  set  is  used.  Control  returns  to  DATGEN. 
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HQPT4 


This  subroutine  computes  the  PD,  QP  and  VTDD  values  from  equations, 
depending  on  the  time  T.  These  equations  are  used  because,  for  this 
scenario,  continuous  input  values  in  distinction  to  discrete  values 
are  pertinent.  After  the  PD  and  QP  values  have  been  determined  in 
degrees/sec2,  they  are  converted  to  radians/sec2.  Control  returns  to 
DATGEN. 

ROLSWA 

This  routine  implements  the  (LAMARS)  drive  logic- (human  pilot)  sensory 
apparatus  simulation  for  the  roll  and  sway  axes.  In  addition,  it  computes 
various  performance  measures,  stores  them  on  file  and  prints  them  on  a 
line  printer.  The  routine  reads  the  data  file  output  from  a scenario 
generation  (QA,  PA,  RA,  UPOD,  VPQD,  WPQD)  and  this  data  is  used  as  input 
to  the  drive  logic-sensory  apparatus  simulation.  The  various  parameters 
(constants)  used  in  the  drive  logic  simulation  are  input  in  the  initiali- 
zation section.  Also,  a number  indicating  the  type  of  drive  logic  (linear, 
subliminal  or  Parrish)  is  read  into  LOGIC,  as  well  as  a device  number  from 
which  the  scenario  data  are  to  be  read  (NDEV).  Three  passes  at  the  drive 
logic-sensory  apparatus  simulation  equations  are  necessary  to  collect  the 
data  for  the  performance  measures.  In  the  first  pass,  the  angular  velocity 
axis  gain  is  set  to  zero.  In  the  second  pass,  the  specific  force  axis  gain 
is  set  to  zero.  The  third  pass  simulates  the  drive  logic  using  the  input 
data  values  for  both  the  angular  velocity  and  specific  force  axis  gains. 

After  the  state  derivatives  (ZD)  are  computed,  the  values  to  be  used  later 
for  the  performance  measures  calculations  are  accumulated  in  the  AXIS 
array.  AXIS  is  dimensioned  to  allow  for  two  axes,  24  measures  and  eight 
different  scenarios  or  parameter  sets.  The  state  derivatives  are  "inte- 
grated" using  a second-order  Adams- Bashforth  numerical  approximation 
(see  Appendix)  to  obtain  the  states  (Z).  This  process  continues  for  each 
pass  until  the  maximum  time  (TMAX)  is  reached.  The  mean  of  each  AXIS  value 
is  computed  and  the  desired  number  of  measures  columns  per  page  is  read. 

This  entire  process  is  repeated  from  the  top  of  the  routine  for  new  parameters 
or  a new  scenario  (or  both)  until  the  desired  number  of  columns  is  reached 
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(NCOL  = LOOP).  The  performance  measures  are  then  computed  (including 
a change  in  units  from  radians  to  degrees  and  application  of  various 
normalizations),  written  to  a file  and  printed  on  the  line  printer. 

If  additional  runs  are  desired,  a 0 is  input;  otherwise  a —1  indicates 
the  end  of  computation. 

PITSUR 

This  routine  simulates  the  (LAMARS)  drive  logic  and  (human  pilot)  t 

sensory  apparatus  for  the  pitch  and  surge  axes,  and  also  computes  per- 
formance measures.  The  program  logic  is  identical  to  that  used  in  the 
routine  ROLSWA,  with  obvious  changes  for  pitch  and  surge  axes  except: 

• No  provision  is  made  for  other  than  linear  drive 
logic . 

YAWHEV 

This  routine  simulates  the  drive  logic  and  sensory  apparatus  for 
the  yaw  and  heave  axes.  It  uses  the  same  logic  used  by  the  routines 
ROLSWA  and  PITSUR  (with  the  obvious  changes  to  account  for  the  [uncoupled] 
yaw  and  heave  dynamics)  except: 

• No  provision  for  nonlinear  drive  logic,  but  yaw 
washout  may  be  first-  or  second-order. 

• Only  one  pass  is  needed  for  drive  logic-sensor 
computations,  since  yaw  and  heave  axes  are 
uncoupled. 

THRES 

This  routine  imposes  the  complementary  function  to  the  indifference 
threshold  which  is  used  in  the  subliminal  (nonlinear)  washout  scheme. 

This  nonlinear  washout  is  used  in  the  sway  axis  only  so  this  routine  is 
called  by  ROLSWA.  At  its  conclusion,  it  returns  control  to  ROLSWA. 

THRESD 

This  routine  is  also  called  by  ROLSWA,  to  impose  the  partial  deriva- 
tive of  the  complementary  function  to  the  indifference  threshold  with 
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respect  to  its  argument.  It,  too,  is  part  of  the  Subliminal  washout 
scheme.  Further  explanation  of  the  functions  computed  by  THRES  and 
THRESD  can  be  found  in  Appendix  B of  the  Engineering  Flight  Simulation 
Validation  final  report. 

This  concludes  the  descriptions  of  the  subroutines  which  comprise 
this  software  package.  Listings  of  each  subroutine  can  be  found  on  the 
pages  which  follow. 
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C 

C MAINLINE  FOR  PLOTTING 

C 

C SETS  UP  VARIABLE  NAMES f SCALE  FACTORS 

C SETS  UP  MAX  TIME  AND  TIME  SCALE  FACTOR 

C CALLS  TIME  HISTORY  GENERATION  SUBROUTINE 

C 

DIMENSION  NVLA( 150) 

COMMON/OUTPUT/NV  t KVLA ( 21 > f SF  < 21 > f KV  < 20 ) 

COMMON/ 10/ 1 IN f IOUT » I TERM  f I PLOT  f I ME AS 1 r IMEAS2 
C 

C INITIALIZE 

C 

DATA  (NVLA(I)f I-1f150)  /150*1H  / 

DATA  <KVLA(I>fI=1f21>f(SF<I>fI=1f21>  /21*1H  f21*0.0/ 

DATA  <KV(I)fI=1f20)fNV  /21*0/ 

DATA  IINfI0UTfITERMfIPL0TfIMEAS1fIMEAS2  /20f21 f5f5f22f23/ 

C 

C INPUT  LIST  OF  POSSIBLE  VARIABLE  NAMES 

C 

DO  20  1=1 f25 

MIN=6*(I-1)+1 

MAX=6*I 

READ(IINfIO)  (NVLA(J)fJ=MINfMAX) 

10  FORMAT  <5(A5f1X)fA5) 

20  CONTINUE 

C 

C READ  NUMBER  OF  VARIABLES f VARIABLE  NAMES f 

C SCALE  FACTORS f MAX  TIMEf  TIME  SCALE  FACTORS 

C 

READ* I IN f 30)  NV 
30  FORMAT (12) 

IFCNV.LE.20)  GO  TO  40 
CALL  ERRMES(l) 

GO  TO  1000 
40  DO  50  1=1 fNV 
50  READ< I IN f 60 ) KVLA< I > fSF< I > 

60  FORMAT (A5f 1XfF10.4) 

READ(IINf70>  TMAXTfTIMSF 
70  FORMAT ( 2F 1 0 • 4 ) 

SF  < NV+ 1 ) =T I MSF 
C 

C CHECK  NAMES 

C 

DO  100  1=1 fNV 
DO  80  J=i f 150 

IF ( KVLA ( I > . EQ.NVIA* J) ) GO  TO  90 
80  CONTINUE 

CALL  ERRMES<2> 

GO  TO  1000 
90  KV< I )=J 

100  CONTINUE 

C 

C RUN  TIME  HISTORIES 

C 

CALL  SENGEN  < TMAXT ) 

C 

C THE  END 

C 

1000  STOP 
,END 
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C*******  **********  ****************************************** 

C*  **************************************  ******************* 

SUBROUTINE  SENGEN < TMAXT ) 

C 

C MAIN  ROUTINE  TO  GENERATE  SCENARIO 

C 

DIMENSION  Z<9>  .ZD<9> ,Z0(9> »ZB0<9> .Y<9> . YD<9) .XAUX (3) 
COMMON/SCENE/  NSCEN 

COMMON/COMDAT/  NDISC . TMAX < 50 ) . P ( 50 ) .0(50 > . VT ( 50 ) 

COMMON/IO/  I IN. IOUT. ITERM. IPLOT. IMEAS1 , IMEAS2 
C 

C INITIALIZE 

C 

DATA  <y<I).I=1.9). <YD<I).I=1.9>,<XAUX<I>,I=1.3)  /21*0.0/ 

DATA  N.IT. IDT. INDEX. DT  /9.0.1.1..1/ 

C 

C READ  SCENARIO  NUMBER  AND  DISCRETE  DATA  VALUES 

C 

READ! I IN . 10)  NSCEN 
10  FORMAT (12) 

IF(NSCEN.EQ.3>  GO  TO  30 
READ (IIN.10)  NDISC 
DO  15  1=1. NDISC 

15  READ( I IN. 20 ) TMAX( I ) , P< I ) .0< I ) . VT < I ) 

20  FORMAT (4F10 . 4 ) 

C 

C CALL  KINEMATICS  SUBROUTINE.  WRITE  DATA 

C 

30  T=IT/10. 

CALL  DATGEN < T . DT . INDEX . PA . QA » RA . UPOD . VPOD . WPOD . Z . ZO . ZD . ZDO ) 

WR I TE ( IOUT . 40 ) P A . OA . R A . UPOD . VPOD . WPOD 
40  FORMAT (6(F11.5»1X>) 

C 

C CHANGE  UNITS.  CALL  PLOTTING  ROUTINE 

C 

DO  50  J-1.6 
Y(J)=Z< J)*57«3 
50  YD< J >=ZD ( J) #57 . 3 

Y(9)=Z<9)*57.3 
YD(9)»ZD<9>*57.3 
DO  60  J=7 . 8 
Y < J ) =Z  < J > 

60  YD ( J)=ZD<  J> 

XAUX  < 1 ) =UP0D/32 « 2 
XAUX  < 2 > =VF'0D/32  • 2 
XAUX  < 3 ) =WP0D/32 . 2 

CALL  OUTR ( T . Y . YD . N . XAUX . TMAXT . IEND . IOCNT ) 

C 

IFIT.GT. TMAXT)  GO  TO  70 
IT=IT+IDT 
GO  TO  30 
70  RETURN 
END 

C*****«*********«*************************»*******************»***** 

c************************************************************* 

SUBROUTINE  DATGEN ( T . DT . INDEX , PA . OA . RA . UPOD . VPOD . WPOD . Z . ZO . ZD . ZDO ) 

C 

C IMPLEMENTS  KINEMATIC  EQUATIONS 
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COMMON/SCENE/  NSCEN 

DIMENSION  Z <9  ) »ZD(9 ) »ZO<  9 > r ZD0<9) r OT<  2) r ISCEN< 3) *ZI < 9) 
DATA  0TU),0T<2)  /2*0.0/ 

G=32.2 

XLX=17.67 

XLZ=-3.1 

IF(T.NE.O.O)  GO  TQ  50 
INITIAL  CONDITIONS 


DO  40  J=l»9 

ZI < J)=0.0 
Z<J>=0.0 
Z0<J)=0.0 
ZD0(J)=0.0 

40 

ZD< J>=0.0 

A0=  < 07 

IF ( NSCEN «GT .4) 

ZI (5)=«07 

IF(NSCEN.E0.4) 

ZI<7)«420. 

IF  < NSCEN . EO  < 4 > 

ZI(8)=-5. 

IF  < NSCEN >LT  «4) 

ZI(7)=560. 

IF  < NSCEN. GT .4) 
DO  45  J=1 r? 

ZI(7)-800. 

45 

Z(J>=ZI(J) 

C 

50 

DELT2=DT/2. 

DO  60  J“1 »9 
ZO(J)»Z<J> 

60 

ZDO(J)=ZD<J) 

C 

C 

GET  DATA 

C 

IF ( NSCEN • NE • 3 ) 

CALL  READAT (Tr INDEX r PD rQDrVTDD) 

IF ( NSCEN • EQ  < 3 ) 

CALL  HQDT4(TrZ'PD'QDrVTDD> 

C 

C 

KINEMATIC  EQUATIONS 

C 

C 

70 

ZD< 1 )=PD 
ZD(2)=QD 

IF  < NSCEN • LE • 2 ) 
ZD(4)-Z(1) 

ZD<2)«G*Z<1)*SIN<Z<4>>*<1./C0S<Z<4>>**2+1>/Z<7> 

ZD<5)-Z<2>*C0S(Z<4))-Z(3)*SIN<Z<4>) 

ZD(6)=Z(2)*SIN<Z<4) >+Z(3)*C0S(Z<4) ) 

ZD(7)=Z<8> 

ZD(8)=VTDD 

IF ( NSCEN. GT .4)  GO  TO  74 
0T<1>«Z<8)+G*Z<S> 

0T(2)=-Z(7>*Z<2>-G»C0S(Z(4>) 

GO  TO  76 

ZD<9>-Z<2>-G*<1.-C0S<Z<4>)>/Z<7>-Z<9)*Z<7)/ZI<7> 

0T< 1 )«Z(8)+G9Z<5)+Z<7)*< A0+Z<9> )*<Z<2>-ZD<9)> 
0T<2)“-Z<7)*Z(2)-G*C0S<Z<4>  >+ZD<9)*Z<7)+Z<8)*<A0+Z<9> ) 
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76  PA=Z ( 1 ) 

QA=Z<2) 

RA=Z  < 3 > 

UPOD=OT  < 1 ) +XLZ*ZD  < 2 ) +XLZ*Z ( 3 ) *Z  < 1 ) -XLX*Z  < 2 > **2-XLX*Z ( 3 ) **2 
VPOD=-XLZ*ZD  < 1 ) +XLX*ZD ( 3 ) +XlZ*Z ( 3 > *Z  < 2 > +XLX*Z  < 1 > *Z  < 2 > 
UPOD=OT < 2 > -XLX*ZD ( 2 > -XLZ*Z  < 1 ) **2-XLZ*Z  < 2 > **2+XLX*Z  < 1 ) #Z  < 3 ) 
C 

C INTEGRATE 

C 

DO  80  J=1 » 9 

80  Z < J ) =Z0  < J ) +DELT2*  <3.*ZD(J)-ZD0(J)> 

Z(3)=G*SIN(Z<4>)/Z<7)+Z(1)*<Z<?)+A0) 

IF < NSCEN . LE . 4 > Z<3)=G*SIN<Z<4) >/Z<7) 

IF ( NSCEN • NE • 3 > GO  TO  90 
IF ( T . GT • 6 . 0 . AND . T . LE .16.5) 

+ Z(1  ) = <C0S<Z<4) ) )**2/(3.'*SIN(Z<4)  >> 

IF(T.GT.5.0.AND.T.LE.17.0) 

+ Z<4>=AC0S(l./(2.+(T-5.0)/3.>) 

C 

90  RETURN 

END 

C**********************************************»******* 
c******»******************************************************** 
SUBROUTINE  READAT < T » I ? PD f QD » VTDD > 

C 

C READ  DISCRETE  SCENARIO  GENERATION  VALUES 

C 

COMMON/COMDAT/  NDISC»TMAX(50)»P(50)rQ(50),VT(50> 

N=NDISC 

C 

IF  < T . GE • TMAX  < I ) > 1 = 1 + 1 
PD=P< I > 

QD=Q< I ) 

VTDD=VT  < I > 

C 

RETURN 

END 

C***<t***********«*************»**************«***********0 
c»*»*»*****************************»*********************»** 
SUBROUTINE  HaDT4<TrZrPDrQD,VTDD> 

DIMENSION  Z<8) 

DPR  = 57.3 
G = 32.2 
VT  = 560.0 
XNZD  = 1./3. 

C 

PD  = 0.0 
QD  = 0.0 
VTDD  = 0.0 
C 

IF < T . GT . 2. 0. AND . T • LE . 3.0 > PD  = 30.0 
IF < T .GT • 4. 0. AND. T • LE . 5.0 ) PD  = -30.0 
IF(T.GT.5.0. AND. T « LE « 6 • 0 ) PD  = 5.5 

IF<T.GT.6.0.AND.T.LE.16.5>  PD  = -DPR*Z ( 1 )*XNZD*COS (Z < 4) ) * 

+ <1.  + <SIN<Z<4>>  >#*2. )/(SIN(Z(4>) >**2. 

IF( T • GT • 16. 5. AND. T .LE .17.5)  PD  = -0.57 
C 

IF<T,GT.2,0, AND. T ,LE .5.0) 


11 


non 


PAGE  1-4  SENLIB  FRX  09-JUN-78  16J57 

+ QD  = DPR*G*Zd)*SIN<Z<4>>*d.+l./<C0S<Z<4>>>#*2.)/VT 
IFd.GT .5.5. AND. T ,LE. 17.0) 

+ QD  = DPR*G*XNZD*d.  + <C0S<Z<4>))**2.)/VT 
C 

PD  = PD/DPR 
QD  = QD/DPR 
RETURN 
END 

C********************************************************** 

c»********************************«************************* 

SUBROUTINE  ERRMESdERR) 

ROUTINE  TO  PRINT  OUT  ERROR  MESSAGES 

COMMON/IO/IINr IOUTf ITERMf IPLOTf IMEAG1 . IMEAS2 
IF<  IERR . EQ . 1 ) WRITEdTERMrlO) 

IF ( IERR.EQ.  2)  WRITEdTERMf  20) 

10  FORMAT <5Xf 'NO  MORE  THAN  20  PLOTTING  VARIABLES  PERMITTED  ') 

20  F0RMAT(5Xf 'PLOTTING  VARIABLE  NAME  NOT  RECOGNIZED  ') 

RETURN 

END 
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C***«OUTR-START 

SUBROUTINE  OUTR<T p Z» ZD » NDIM r XAUX* TMAXT p IENDp IOCNT) 

DIMENSION  Z( 1 ) f ZD< 1 ) pXAUX< 1 > pXOUT ( lOOf 21 ) 

DIMENSION  IDATE(2) 

C 

C MODIFIED  FOR  GPDPLT  GLT  19-SEP-1977 

C 

COMMON/OUTPUT/  NV p K VL A < 2 1 ) p SF  < 2 1 ) p K V < 20  > 

C 

C 

C 

IF<T.NE.O. >00  TO  30 
IOCNT-O 
I0MAX-100 
NUT-NV+1 

C 

C PLOT  AXIS 

C 

CALL  GPDPLT  < XOUT  p NU , I OCNT  r KVLA t SF » TMAXT » IEND ) 

C 

C*#**********************************»******************************** 

c 

C UPDATE  COUNTER 

C 

30  I0CNT=I0CNT+1 
C 

C STORE  DATA 

C 

XOUT  < IOCNT  »NVT>=T 
C 
C 

DO  35  J=1,NV 
I=KV(J) 

IF(I.GT.50)  GO  TO  32 
XOUT  < IOCNT  » J)*Z< I > 

GO  TO  35 

32  IFd.GT. 100)00  TO  34 

XOUT ( IOCNT  » J ) =ZD ( I -50 ) 

GO  TO  35 

34  XOUT  < IOCNT » J)=XAUX( 1-100) 

35  CONTINUE 
C 

C TIME  TO  PLOT  ? 

C 

IFd.GT . TMAXT  > IEND=1 

IF(IEND.NE.O) GO  TO  40 

IF( IOCNT .LT . IOMAX)GO  TO  100 

c 

C PLOT  SEGMENT 

C 

40  CALL  GPDPLT<XOUTfNUpIOCNT»KVLA*SF»TMAXT»IEND) 

C 

C UPDATE  OUTPUT  MATRIX  AND  RESET  COUNTER 

C 

DO  50  J=1»NVT 

50  X0UT(1»J>»X0UT<I0CNT»J) 

I OCNT- 1 

c 


I 
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100  RETURN 
END 

C*i**OUTR-END 

SUBROUTINE  GPDPLT ( DATA  ? NVAR?  NPTS ? LABELS ? SF ? TMAXT ? IEND) 

DIMENSION  DATA( 100 r 21 ) r LABELS (21)?SF(21)?IDATE(2) 

C 

C GENERAL  PURPOSE  DATA  PLOTTING  GLT  3/1/77 

C MOD  FOR  XV15  GLT  AUG/77 

C 

C DATA  AN  ARRAY  OF  NPTS  OF  NVAR  VARIABLES  TO  BE  PLOTTED  VERSUS 

C THE  INDEPENDENT  VARIABLE?  DATA ( I ? NVAR+ 1 ) ? I=1?NPTS. 

C LABELS  A CORRESPONDING  SET  OF  5 CHARACTER  VARIABLE  NAMES. 

C SF  A CORRESPONDING  SET  OF  SCALE  FACTORS  (UNITS/INCH). 

C TMAXT  THE  ULTIMATE  MAXIMUM  VALUE  OF  THE  INDEPENDENT  VARIABLE? 

C WHICH  IS  ONLY  USED  TO  DETERMINE  THE  PAGE  WIDTH. 

C IEND  = 0 MORE  DATA  TO  COME 

C 1 THE  LAST  SET  OF  DATA 

C 

C GPDPLT  WILL  PLOT  THE  DATA  USING  A ZETA  PLOTTER!  5 VARIABLES  TO 

C A PAGE.  A PAGE  IS  11  INCHES  HIGH  AND  A MULTIPLE  OF  8.5  INCHES 

C WIDE.  THE  WIDTH  IS  DETERMINED  BY  TMAXT  AND  SF<NVAR+1>.  THE 

C INITIAL  PEN  POSITION  SHOULD  BE  CENTERED  (5.5  INCHES  FROM  THE  TOP) 

C AND  1 INCH  TO  THE  RIGHT  OF  THE  LEFT  EDGE. 

C 

C GPDPLT  IS  STRUCTURED  TO  BE  CALLED  SEVERAL  TIMES  TO  COMPLETE  A 

C SET  OF  PLOTS.  ON  THE  FIRST  CALL  (NPTS  = 0)?  GPDPLT  IS  INITIALIZED 

C AND  AXES  ARE  PLOTTED.  ON  SUCCEEDING  CALLS  (NPTS.NE.O)  WHATEVER 

C DATA  ARE  AVAILABLE  (1  = 1, NPTS)  ARE  ADDED  TO  THE  PLOTS.  ON  THE 

C FINAL  CALL  (IEND  = 1)  THE  LAST  SET  OF  DATA  IS  PLOTTED!  A 

C LABEL  IS  APPENDED  TO  THE  END  OF  EACH  VARIABLE  PLOT!  AND  THE 

C PAPER  IS  EJECTED. 

C 

COMMON/ 10/ 1 IN? IOUT  ? I TERM  ? IPLOT  ? IMEAS1 ? IMEAS2 
DIMENSION  PDATA( 1024) 

C 

C*********************************************************************** 

c 

IF  ( NPTS . NE . 0 ) GO  TO  200 
C 
C 

C *********************************************************************** 

c 

C INITIALIZE 

C 

NTV-NVAR+1 
XIN=DATA( 1 ?NTV> 

XISF=1.0/SF(NTV> 

XMAX=TMAXT*XISF 

PWIDTH  = 8.5  * ( 1 + INT (XMAX/8.5)  ) 

IF (PWIDTH-XMAX.LT • 1 .0)PWIDTH=PWIDTH+8.5 
IDIR  = -1 


C 

C ************************************************************ ******  *** 

c 

C PLOT  AXES 

C 


CALL  PLOTS (PDATA? 1024? IPLOT) 
DO  50  I • 1 ? NVAR 
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IVAR=I 

IF<I.GT.5.AND.I.LT.11)IVAR  = MINO< 14-1 fNVAR+6-I ) 
IF(I.GT.15.AND.I.LT.21)IVAR  = MIN0(36-I fNVAR+16-I ) 

IM1  = I VAR  - 1 
M0D5=  MOD  ( IM1 r 5 ) 

YORG  = 3.5  - 1.5  * MOD5 
XORG  = < IM1/5)  * PWIDTH 
C 

SCALE  ■ 0.5*  SF(IVAR) 

LFRAC  = 2 - INT ( ALCG 1 0 ( ABS ( SCALE  > ) > 

IF (LFRAC.LT .0)LFPAC  - 0 
C 

CALL  SYMBOL ( XORG-0 . 9 1 YORG t « 14 r LABELS < I VAR) r 0. r 5) 

C 

50  CONTINUE 
GO  TO  450 
C 
C 

C *********************************************************************** 

c 

C PLOT  DATA 

C 

C PLOT  PAGES  LEFT  TO  RIGHT  IF  IDIR  = +1 

C PLOT  PAGES  RIGHT  TO  LEFT  IF  IDIR  = -1 

C 

200  DO  400  II=1fNVAR 
I=II 

IF< IDIR.EQ.-l ) I=NVAR+1-II 
I VAR=I 

IF  < I . GT . 5 . AND . I . LT . 1 1 ) I VAR=MINO  < 16-1 f NVAR+6-I > 

IF  < I . GT . 15. AND .I.LT.21 ) I VAR=MINO<  36-1 fNVAR+16-I ) 

IM1=IVAR-1 
M0D5CM0D< IM1 f5) 

Y0RG=3.5-1.5*M0D5 
X0RG=<IM1/5)#PUIDTH 
LDIR=(  1 - 2 * M0D(M0D5f2>  > * IDIR 
C 
C 
C 

c 
c 
c 
c 
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PLOT  THE  IVAR'TH  VARIABLE?  FOR  LDIR  = 1.  LEFT  TO  RIGHT 

-If  RIGHT  TO  LEFT 

YISF=1 ./SF< I VAR) 

LO  = 1 

IF (LDIR.EQ.-l )LO=NPTS 
X=XORG+ < DATA  < LO f NTV > -XIN ) *XISF 
Y=YORG+DATA (LOf I VAR  > *YISF 
IF(Y.GT.5.0)Y=5.0 
IFIY.LT.-5. >Y=-5.0 

IF< IEND.NE.O. AND. LDIR.EQ.-l )CALL  SYMBOL <X+ . 1 fY- ,05f . 14 t 
+ LABELS! I VAR) fO. f5) 

CALL  PLOT (Xf Yf3) 

DO  360  L =2fNPTS 
LO=LO+LDIR 

X=<DATA(LOfNTV)-XIN>*XISF 
Y= YORG+DATA < LO  f 1 VAR ) * YI SF 
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IF(Y.GT.5.0)Y«5.0 
IF(Y.LT.-5. >Y=- 5.0 
CALL  PLOT (X+XORGf  Yf2> 

360  CONTINUE 

IF(  IEND.NE.O.  AND . LD  IR . EO . 1 ) CALL  SYMBOL  < X+XORG+ . 1 f Y- . 05 1 

+ • 14 f LABELS ( I OAR) fO. r5> 

C 

C 

400  CONTINUE 
C 
C 

IDIR=-IDIR 

C 

IF< IEND.EQ.0 )GQ  TO  450 
DO  420  I =1 fNVAR 
I VAR=I 

IF<  I . GT . 5.  AND.  I .LT  .11)  I VAR  - MIN0U6-I  fNVAR+6-I  > 

IF( I . GT . 15 . AND .I.LT.21) I WAR  - MIN0<36-I fNVAR+16-I > 

IM1  - I VAR  - 1 
M0D5=  MOD< IM1 .5) 

YORG  = 3.5  - 1.5  * MODS 
XORG  - ( IM1/5)  * PWIDTH 
C 

SCALE  » 0.5*  SF(IVAR) 

LFRAC  = 2 - INT ( AL0G1 0 < ADS ( SCALE ) > ) 

IF  < LFRAC . LT . 0 > LFRAC  - 0 
C 

CALL  NUMBER  < XORG-0 . 8 f YORG+0 . 5 r . 1 4 f SCALE f 0 . r LFRAC  ) 

CALL  PLOT (XORGrYORG+O.5'3) 

CALL  PL0T<X0RGfY0RGf2> 

CALL  PLOT ( XORG-O • 1 » YORG f 3 ) 

CALL  PLOT < XORG fY0RGf2> 

CALL  PLOT <X0RGfY0RG-0.5f2> 

C 

IF < I VAR . NE . NVAR . AND . YORG . NE . -2 . 5 ) GO  TO  420 
C 

CALL  AXIS ( XORG r -4 . 1 4HTIME  f -4  f 1 . fO.  fXINfSF(NTV)  ) 

420  CONTINUE 

CALL  PLOT < < NVAR/5 ) *PWIDTH+8 . 5 fO f 3 ) 

^50  CALL  PLOTE 
500  RETURN 
END 
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DIMENSION  Z<18)fZD(18)f ZOO ( 18 ) fZO( 18 ) fXK ( 2) f SCALE < 2 > 

DIMENSION  AXIS  < 2 f 29*8) f XMEAS  < 8 f 24 ) f ILABEL (24)fUNIT(2) 

DIMENSION  LUN(2>  fXAUX<8)  f Y ( 18  > f Y1.K18 ) r PEAK ( 2 ) fSCALE2 ( 2 ) fXK2<2> 
DATA  IIN.ITERMf1MEAS1fIMEAS2  /20f5f21f22/ 

DATA  < ILABEL ( I)fI-1f24)  /5HPDLA  fSHPDLI  f5HPDL2  f5HPDL3  f 
+ 5HPDL4  fSHPDLS  f 5HPSENA  f 5HPSEN1  f 5HPSEN2  f5HPSEN3  f 5HPSEN4  f 

+ 5HPSENS 1 5HPNL A fSHPNLI  fSHPNL2  f5HPNL3  f5HPNL4  fSHPNLS  f 

+ SHGS1  f 5HGS2  fSHGSS  fSHPACC  fSHPVEL  f5HPDIS  / 

LOOP=l 

r 

U 

C GET  SCENARIO  FILENAME 

l: 

10  READ < I IN f 20)  NDEV 

20  FORMAT  < 12) 

C 

C CONSTANTS 

C 

READ  < I IN  f 30 ) TMAXfDT 
30  FORMAT ( 2F9  * 4 ) 

DELT2=DT/2. 

READ  ( I IN  F 70 ) GP  F WP  F U1R  F W2R  f W3R  f U4R  f G2R  f C2  Y 
70  FORMAT  ( 8F9 • 4 > 

READ ( I IN f 76 ) LOGIC f BP f BK fC2YP 
76  FORMAT (I2f3F9*4) 

C 

C 3 PASSES  AT  DATA 

C 

DO  160  IPASS=1f3 
XKP-GP 
XK2R=G2R 
XKY=1.0 

I F < I PASS  f EQ « 1 ) XKP=0.0 
IF  < I PASS f EG • 2 ) XK2R=0.0 

C 

C INITIALIZE 

C 

DO  80  1=1 » 18 
Z(I)=0.0 
80  ZD<I)=0.0 

Z(18)=XKP 
IT=0 
IDT=1 
t: 

C DATA  COMPUTATION  LOOP 

C 

90  T=IT/10. 

DO  100  1=1 f 18 
ZO(I)=Z<I) 

100  ZDO ( I ) =ZD ( I > 

C 

READ < NDEV rllO)  PAfOA.RAf UPOD f VPOD f WPOD 
110  FORMAT (6<F11f5f1X)) 

VP0D=VP0D/32.2 

IF< LOGIC . FQ f 2)  GO  TO  112 

C 

ZD  < 1 ) =XKP*PA-WP*Z  < 1 ) 

ZD<  2 >=ZD  < 1 ) -Z<  4 ) 

GO  TO  114 


! 

i 
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.112  Z ( 1 8 ) “ Art  I N 1 ( Z ( 1 8 ) .XKP>  I AMAX1 (Z( 18) rO.  >-Z<18> 

ZD  < 1 > -Z ( 1 8 > #F'A-UP#Z  < 1 ) 

ZIU2)=ZD(1)-Z(4> 

ZD( 1 7 ) =PA~UP*Z  ( 17  > 

ZD( 18 )= ( PA-ZD ( 1 ) >*ZD(17)-BP#Z(1 >*Z(17) -BK*(Z( 181-XKP) 

ZD  v 18  > -AMAX1 < ZD< 18)f-0.5> 

114  continue; 

ZD  < 3 > “XN2R*VP0D+Z < 2 > 

ZD(4)=U4R**2*ZD(3)+WlR*U4R#*2*Z(3)-< 1 . 4*W3R+W2R >*Z < 4 ) 
v IF  ( LOGIC* NE  • 3 > ZD ( 5) =ZD ( 3 ) -1 . 4#C2Y#» . 5*Z < 5 > ~C2Y*Z ( 6 > 

\ IF(L08IC.EQ.3)  ZD (5 >=ZD< 3) -THRES < . 1 f 1 . 4*C2YP#*. 5*Z (5 ) 

+ +C2YF’*Z  < 6 ) ) *T  HRE3D  (.If  XK2RKUP0D > -1 . 4*C2 Y** . 5*Z  < 5 ) ~C2 Y*Z  < 6 ) 
ZD(6)-Z(5> 

PAF-ZD  (2) 

VP0FD~ZD(5>-Z(2> 

ZD(7)=PA-.2*Z<7) 

ZD(9)=PAF-.2*Z<9> 

ZD( 1 1 ) -1 . 5KUP0D+ . 12*Z  < 1 2 ) - 1 . 7*Z (11)-. 3*Z ( 13 ) 

ZD<12)=WPOD 
ZD( 13)=Z( 11 ) 

ZD<i.4)«1.5*UP0FD+.12*Z<15)-1.7*Z(14>-.3*Z(16> 

ZD< 15>=VP0FD 
ZD( 16)=Z( 14) 

C 

PE=XKP*PA-PAF 

PPE«XKP#Z(8)-Z<10> 

F'AI“Z  ( 8 > -THRES  ( . 035  f Z ( 8 ) > 

PAFI‘-Z ( 10 ) -THRES ( • 035 r Z( 10 ) ) 

PIE~XKPX<PAI-PAFI 
UED=XKY# ( Z(  2 ) +VPOD  > -ZD (5 ) 

VED0=XK2R*VP0D-UP0FD 
VPED-XK2R#Z < 1 1 > -Z  < 14  > 

VP0ID«Z(in-THRES<0.1fZ(ll>) 

VPOFI D=Z (14) -THRES ( 0 . 1 . Z < 14 ) ) 

UIED--XK2R*VP0ID-VP0FID 

C 

C PASS  1 

C 

IF(IPASS.NE.l)  GO  TO  120 

AXIS(2f2fL00P)=AXIS(2f2»L00P)+UED**2 

AXIS (2  f 4 f LOOP ) “AXIS ( 2 r 4 f LOOP ) + (Z( 2 ) +UP0D  > #ZD ( 5 ) 

AXIS ( 2 r 8 f LOOP  > = AXIS ( 2»8r LOOP ) +UPED##2 
AXIS ( 2 f 1 Of LOOP  > = AXIS  < 2 r 10  f LOOP )+Z(ll)*Z(14) 

AXIS(2f 13f LOOP )»AXIS(2f 13f LOOP >+VIED**2 
AXIS(2fl5f LOOP )“AXIS(2f 15 f LOOP ) +UPOID#UPOFID 
AXIS ( 2 f 25  f LOOP ) =AXIS ( 2 • 25  r LOOP ) +VPOD*MPOFD 
AXI5<  2 f 28  f LOOP ) =AXIS ( 2 r 28f LOOP )+VED0##2 
C 

C PASS  2 

C 

120  CONTINUE 

IF( IPASS .NE . 2 ) GO  TO  130 
AXIS(lf2f LOOP )=AXIS(l»2f LOOP >+PE**2 
AXIS ( 1 r 4 f LOOP  > “AXIS! 1 »4. LOOP ) +PA*PAF 
AX  I G ( 1 f 8 f LOOP  > =*  AX  I S < 1 . 8 r LOOP  ) +PPE**2 
AXlS<lflOf  LOOP ) =AXIS ( 1 » lOfLOOP  >+Z(8)*Z(10) 

AXIS( 1 r 13f LOOP)=AXIS( 1 r 13rL00P)+PIE#*2 
AXI3( If 15fL00P)=AXIS(lf 15fL00P)+PAI*PAFI 
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C PASS  3 

w 

130  CONTINUE 

IF  ( I PASS  < NE  < 3 ) GO  TO  140 

AX IS  < 1 , 1 f LOOP ) -AX I S < 1 f 1 r LOOP ) +PA**2 

AXIS<1f3.L00P)=AXIS<1f3fL00P)+PE**2 

AXIS ( 1 < 5 f LOOP ) -AXIS (1 , 5 r LOOP > +PA*PAF 

AXIS<1f6fLOOP)=AXIS<1f6fLOOP)+Z<8)#*2 

AXIS<1f7fL.OOP)=AXIS<1f7fLOOP)+Z<10)**2 

AXIS<1f9fL00P)=AXIS<1f9fL00P)+PPE**2 

AXIS<1f11fL00P)=AXIS<1f11fL00F)+Z<8>*Z<10> 

AXIS< 1f12fL00P)=AXIS<1f12fL00P)+PAI**2 
AXIS<1f14fL00P)=AXIS<1f14fL00P)+PIE**2 
AXI S ( 1 f 16  r LOOP ) —AX IS<1f16f LOOP ) +PAI*PAFI 
AXIS<1f21fL00P)=AMAX1 <AXIS<1<21 fLOOP) fABS<ZD<2)*57.3) ) 
AXIS<1f22fL00P)=AMAX1<AXIS<1f22fL00P) f ABS < Z < 2) *57. 3 ) ) 
AX I S < 1 f 23  f LOOP ) “AXIS  < 1 < 23  f LOOP ) +PAF**2 
AX I S < 1 < 24  < LOOP ) “AXIS  < 1 f 24  fLOOP ) +F'AFI**2 
C 

AXIS(2f1f  LOOP  > = AXIS  < 2 f 1 f LOOP > + < Z < 2 ) +UPOD > **2 

AXIS  ( 2 f 3 fLOOP  > =AXIS  (2f  3 fLOOP  ) +UED**2 

AXIS ( 2 f 5 fLOOP ) “AXIS  < 2 f5f  LOOP ) + <Z  < 2) +VPOD ) *ZD  < 5 ) 

AXIS  ( 2f  6 fLOOP  ) “AXIS  ( 2 f 6 f LOOP  > + Z < 1 1 ) **2 
AXIS  < 2f  7fL00P ) =AXIS  < 2 f 7 f LOOP ) +Z  < 14 ) **2 
AXIS  ( 2 f 9fL00F'  ) “AXIS  < 2fV  fLOOP  ) +VPED**2 
AXIS  1 2 f 1 1 fLOOP  ) “AXIS < 2 f 1 1 f LOOP ) +Z ( 1 1 ) #Z (1 4 ) 

AXIS ( 2 f 12f LOOP' ) “AXIS < 2 f 12 fLOOP ) +VP0ID##2 

AXIS ( 2 f 1 4f  LOOP  > “AXIS  < 2f 14  fLOOP ) +VIED**2 

AXIS  C 2 f 16  fLOOP  )“AXIS<2f  16  fLOOP  ) +UPOID#UF'OFID 

AXIS C 2 f20f LOOP) “AMAXK AXIS! 2 f20 fLOOP) f ABS<ZD<5)*32.2> ) 

AXIS ( 2 f 21  f LOOP ) “AMAX1  (AXIS(2f21  fLOOP ) f ABS < Z <5 ) *32  < 2 ) ) 

AXIS  < 2f  22  f LOOP' ) “AMAX1  < AXIS  < 2f 22 fLOOP ) f ABS (Z(6  > #32 . 2 ) ) 

AXIS<2<23<  LOOP)  “AXIS  < 2 f23fL00P)+W0FD**2 

AXIS ( 2 f 24  fLOOP) “AXIS  < 2 f24  f LOOP ) +VP0FID**2 

AXIS(2f26f  LOOP' ) “AXIS  < 2 f 26  fLOOP )+VP0D**2 

AX IS  < 2 f 27  f LOOP ) “AXIS < 2 f 27  f LOOP  > +UPOD*VPOF D 

AXIS  < 2 f 29  fLOOP  >“AXIS<2f  29fL0GP  >+UED0**2 

r. 

i.  INTEGRATE 

p 

140  DO  ISO  1=1 f 18 

Z < I ) “ZO ( I ) +DELT2*  < 3 . *ZD (I ) -ZDO  < I ) ) 

IF<I<EQ<8<0R<1< EO • 10 ) Z < I ) “ . 37#Z0< I ) + . 63*ZD < 1-1 > 

150  CONTINUE 
IT=IT+IDT 

IF < T . LT < TMAX  > GO  TO  90 
160  CONTINUE 

i; 

L COMPUTE  MEANS 

C 

N=TMAX/DT 
DO  170  1=1 f2 

DO  170  J=1f29 

IF<J<GE<20< AND < J< LE < 22)  GO  TO  170 
AXIS<IfJfL00P)=AXIS(IfJfL00P)/<N+1> 

170  CONTINUE 


I 


19 


ROLSUA 
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I'AUE  1-4 

IF (LOOP.NE. 1 ) GO  TO  200 
READ (I IN t 190 > NCOL 
190  FORMAT (II) 

200  L00F=L00P+1 

IT (LOOP. LE. NCOL)  GO  TO  10 
C 

C PRINT  MEASURES 


XK  < 1 ) =XKP 
XK  ( 2 ) -XK2R 
SCALE  ( 1 ) =XKP 
SCALE ( 2 ) -XK2R 
UNIT< 1 )=57 • 3 
/ UNIT  ( 2 ) =1 ♦ 0 
‘ LUN ( 1 ) -IMEAS1 
LUN ( 2) =IMEAS2 
PEAK  < 1 ) =XKP 
PEAK  ( 2 )=XKP 
SCALE2 ( 1 ) =XKP 
SCALE2(2)=XKY 
XK2( 1 )=XKP 
XK2<  2 > =XK  Y 
DO  228  1=1 f8 
DO  228  J=1 f24 
228  XMEAS (IfJ)=0.0 

r 

DO  270  IAX=1f2 

IF< SCALE < I AX ) .EQ.O.O)  SCALE ( I AX ) =1 . 0 
DO  230  L00P=1 r NCOL 

XMEAS < LOOP  r 1 ) =AX IS  < I AX  r 1 r LOOP ) *UNI T (I  AX  > **2 

XMEAS  C LOOP  r 2 ) =AXIS( IAXr  2 fLOOP ) /AXIS ( IAX  r 1 r LOOP )/SCALE2 (IAX)##2 
XMEAS(L00Pr4)=(XK2(IAX)*AXIS(IAXf 1 fLOOP >-AXIS( IAX.4rL00P> > 

+ /AXIS  < I AX  r 1 r LOOP )/SCALE2( I AX ) 

XMEAS (LOOP  r3)=(AXIS(IAX » 3 »LOOP ) -AXIS  < IAX  r 2 r LOOP ) ) 

+ /AXIS ( IAX  r 1 fLOOP ) /SCALE2  < I AX ) **2 

XMEAS(L00Pr5)=(AXIS(IAXr 4 »LOOP) -AXIS' IAX»5 fLOOP)  > 

+ /AXISdAXrlr  LOOP ) /SCALE2 ( I AX ) 

XMEAS  < l OOP r 6 > =AXI S ( I AX  r 23  r LOOP ) *UNIT ( I AX  > **2 
XMEAS (LOOP  f7)=AXIS( IAXr  6 fLOOP) 4UNIT (I AX ) #42 
XME  AS  (LOOP r 1 2 ) =AXI  S ( I AX , 7 r LOOP ) 4UNIT  ( I AX  > **2 

XMEAS(LOOP»8)=AXIS(IAXf8fLOOP)/AXIS(IAX»6»LOOP)/SCALE(IAX>**2 
XMEAS ( LOOP  r 10 ) = ( XK ( I AX ) 4AXIS ( IAXr6  r LOOP ) -AXIS ( IAXr 10? LOOP  > ) 

+ /AXIS(  I AX » 6 rl.OOP ) /SCALE  (IAX) 

XMEAS ( LOOP  r 9 > = ( AXIS ( I AX  r 9 r LOOP ) -AXIS ( IAX  r8  r LOOP ) ) 

+ /AXIS(IAXr6rL00P>/SCALE(IAX>**2 

XMEAS ( LOOP  r 11 ) = (AXIS( IAXf 10  r LOOP ) -AXIS ( I AX  r 1 1 r LOOP ) ) 

+ /AXIS ( IAX  * 6 »LOOP ) /SCALE (I  AX ) 

XMEAS  ( LOOP  rl3)=AXIS(IAXr  12  r LOOF* ) 4UNI T (IAX)##2 

XMEAS ( LOOP f 14 >~AXIS( IAXr 13r LOOP >/AXIS( IAXf 12 r LOOP > /SCALE (IAX > **2 
XMEAS ( LOOP  »16)  = (XK(IAX>  4AXIS ( IAXf 12  r LOOP > -AXIS ( IAXr 15  fLOOP ) ) 

+ /AXIS  ( IAXr  12  fLOOP ) /SCALE  ( I AX  ) 

XMEA5(L00P r 15)  = (AXIS (IAXr 14 r LOOP) -AXIS (IAXf 13 r LOOP > > 

+ /AXISdAXf  12  fLOOP) /SCALE  (IAX)  **2 

XMEAS ( LOOP r 17)=(AXIS( I AX r 15 r LOOP ) -AXIS ( I AX r 16r LOOP ) ) 

+ /AXISdAXf  12fL00P)/SCALE(IAX) 

XMEAS  ( LOOP  r 1 8 ) = AX  IS  < I AX  r 24  r LOOP  ) *UN1  T ( I AX  > **2 
XMEAS ( LOUP r 19 >=0.0 
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XMEAS  < LOOP f 20 ) =0  • 0 
XMEAS  < LOOP  f 21 ) =0 • 0 

XMEAS  ! LOOP  f 22)  “AXIS  ! I AX  f 20fLOOF‘  >/F'EAK  ! I AX ) 
XMEAS!LOOPf23)=AXIS!IAXf21fLOOP>/PEAK!IAX> 

XMEAS  C LOOP f 24 ) =AXI  S ! I AX  f 22  f LOOP ) /PEAK < I AX ) 

I F ( I AX . EO . 1 > XME AS  < LOOP  f 22  > =0 . 0 
230  CONTINUE 
C 

DO  260  J=1f24 

WRITE!LUN! IAX) f240>  (XMEAS<NLOOP, J) fNL00F=1 fNCOL) 

240  FORMAT <8F9. 2) 

DO  245  NL00P=1 fNCOL 
TEST =ABS  < XMEAS ! NL00P  » J ) ) 

IF<TEST. GE. 100.0. OR. (TEST. LE. 0.00005. AND. TEST. NE. 0.0))  GO  TO  25 
245  CONTINUE 

UR I TE  < I TERM  f 250 ) I L ABEL ! J ) f ! XMEAS ( NLOOP  fJ) f NLOOP= 1 fNCOL) 

250  FORMAT (1XfA5f1Xf  8F9  *5) 

GO  TO  260 

252  WRITE ( I TERM  1 254 ) ILABEL < J ) > ! XMEAS (NLOOP fJ> fNL00P=1 fNCOL ) 

254  FORMAT  < IX  f AS  f IX  f8E9  • 2 ) 

260  CONTINUE 

WRITE! ITERMi 265) 

265  FORMAT!///) 

270  CONTINUE 

DO  275  L00P=1  fNCOL 

XMEAS ! LOOPf 1 )=AXIS! 2 f 26  f LOOP ) #UNIT ! 2 ) ##2 

XMEAS ! LOOP  f 2 ) = AX IS ! 2 f 28 1 LOOP ) / AX I S ! 2 f 26  > LOOP ) / SCALE  < 2 ) **2 
XMEAS ! LOOP f3)  = !AXIS!  2 f 29 f LOOP ) - AXIS! 2 f28  f LOOP ) ) / AXIS ! 2 f 26 f LOOP ) 
+ /SCALE ! 2 ) #*2 

XMEAS!LOOPf4)=!XK!2)*AXIS!2f26fLOOP)-AXIS!2f25fLOOP) )/ 

+ AXIS ! 2f  26  fLOOP ) /SCALE ! 2 ) 

275  XMEAS! LOOP’ f5)-!AX1S ! 2f 25f LOOF' )— AXIS !2 f27 f LOOP ) >/AXIS!2f26f LQOF* ) 
+ /SCALE <2 ) 

DO  279  J-l f5 

URI TE ! LUN ! 2 ) F 240 ) ! XMEAS ( NLOOP  f J ) f NLOOP= 1 f NCOL ) 

DO  276  NLOOF-1 fNCOL 
TEST=ABS ! XMEAS ! NLOOP  f J ) ) 

IF ( TEST . GE . 1 00 . 0 . OR . ! TEST . LE . 0 . 00005 . AND . TEST . NE . 0 . 

276  CONTINUE 

WRITE! ITERMf250)  ILABEL ! J > f ! XMEAS INLOOPf J ) fNL00P=1 f 
GO  TO  279 

277  WRITE ! I TERM  f 254 ) ILABEL !J > , ! XMEAS ! NLOOP t J ) f NL00P=1 f 
279  CONTINUE 
C 

READ! IINf290)  IANS 
290  FORMAT! 13) 

IF! IANS • EO . -1 ) GO  TO  310 

C 

DO  300  1=1 f2 
DO  300  J=1 f 24 
DO  300  K=1 f 8 
300  AXIS! IfJfK)=0*0 
L00P=1 
GO  TO  10 
C 

310  RETURN 
END 


))  GO  TO  277 
NCOL) 

NCOL) 
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DIMENSION  Z ( 14 ) *ZD ( 14 ) r ZDO <14)*Z0<14)»XK(2)f SCALE ( 2 ) 
DIMENSION  AXIS(2»24f8> fXHEAS (8  f 24) f I LABEL ( 24 > fUNIT < 2 ) 
DIMENSION  LUN (2 ) f PEAK ( 2 ) 

DA  r A I IN  f ITERM  f tilEASl  t IMEAS2  /20f5f21f22/ 

DATA  < I EASEL <I)fI=1f24)  /5MPDLA  f 5HPDL1  f5KPDL2  .5HPDL3  f 
+ 3HP0L4  fSHPDLS  f 5HPSENA f 5HPSEN1 f5HPSEN2f5HPSEN3f5HPSEN4f 
+ 5HP8ENS f 5HPNL A •5HPNL1  f5HPNL2  f5HPNL3  f5HPNL4  f 5HPNLS  f 

+ 5H0S1  f5HGS2  > 5HGSS  f5HPACC  f 5HPVEL  fSHPDIS  / 

LOOF'-l 

C 

i:  GET  SCENARIO  FILENAME  AND  OPEN  FILE 

C 

10  READ ( I IN r 20 ) NDEV 

20  FORMAT (12) 

C 

C CONSTANTS 

C 

READ < I IN f 30)  TMAXfDT 
30  FORMAT  < 2T9 . 4 ) 

DELT2=DT/2. 

READ  < I IN  f 70  > G2PfU0fW3PfU4P 
70  FORMAT ( 4F9 . 4 ) 

C 

C 3 PASSES  AT  DATA 

C 

DO  160  IPASS=1f3 

XK0=0.5 

XK2P=G2P 

IF(IPASS.EQ.l)  XKQ=0.0 
IF ( IPASS.EQ.2)  XK2P-0 • 0 

C 

C INITIALIZE 

C 

DO  80  I=1f14 
Z<I>=0.0 
80  ZD(I)=0.0 

IT-0 
IDT=1 
C 

C DATA  COMPUTATION  LOOP 

t 

VO  T=IT/10. 

DO  100  1=1 f 14 

zo<i>=z<n 

100  ZDO( I >=ZD( I ) 

r 

READ ( NDEV .110)  P A , 0 A r RA f UPOD t VPOD r UPOD 
110  FORMAT (6<F11f5f1X)> 

UP0D=UP0D/32.2 

C 

ZD ( 1 > =XK0*QA-UQ*Z ( 1 > 

ZD(2) *ZD< 1 >+Z<4> 

ZD ( 3 ) =XK2P*UP0D-Z  < 2 ) 

ZD(4)=<U4P**2>*ZD<3>-1.4*U3P*Z<4> 

QAF=ZD<2) 

UP0FD=Z<2> 

ZD(S)=QA-.2*Z(5> 

ZD<7)=QAF-.2*Z<7) 
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y 

1 


ZD<9>=1.5*UP0H+.12*ZC10)-1.7*Z<9)-.3*Z<11> 

ZIK 10)=UP0D 
ZII<11)=Z<9> 

ZrK12)-1.5*UP0FD+.12*Z<13>-1.7#Z<12>-.3*Z<14) 

ZD< 13>=UPQFD 
ZIK14)=Z<i2> 

C 

GE-XNQKQA-QAF 
QPE-XKG#Z  < 6 ) -Z  < 8 > 

GAI-Z  < 6 ) -THRES  < . 035  rZ ( 6 ) ) 

QAFI-Z<  8) -THRES < .035fZ<8>) 

QIE=XKQ*OAI-QAFI 
UEH=XK2P#UPOD-UPOF£l 
UPEU=XK2P*Z  < 9 ) -Z  < 12 ) 

UPOID-Z < 9 ) -THRES  <0.1»Z<9)> 

UPOFID-Z (12) -THRES < 0 . 1 , Z f 12  > > 

UIED=XK2P*UP0ID-UP0FID 

C 

C PASS  1 

C 

IF < IPASS.NE. 1 ) GO  TO  120 

AXIS<2,2,L00P)=AXIS<2»2,L00P>+UEIi**2 

AXIS<2.4,LOOP)=AXIS<2,4»LOOP)+UPOD#UPOFI) 

AXIS  < 2.8 .LOOP) =AXIS<  2, 8 , LOOP  >+UPED##2 
AXIS ( 2 r 10 , LOOP) -AXIS  <2. 10 .LOOP  )+ZC9)#Z<12) 

AXIS  < 2 , 13 .LOOP  > =AXIS<2» 13 » LOOP ) +UIEI)##2 
AXIS<  2 , 15 » LOOP  > “AXIS ( 2 , 15, LOOP  > +UPOID#UPOFID 

C 

C PASS  2 

C 

120  CONTINUE 

IF<IPASS.NE.2>  GO  TO  130 

AXIS ( 1 , 2 » LOOP ) -AXIS  < 1,2, LOOP  > +QE*#2 

AX I S ( 1 , 4 , LOOP ) = AX I S ( 1 » 4 , LOOP ) +QA4QAF 

AXIS  < 1 , 8 f LOOP ) -AXIS  < 1 * 8 r LOOP ) +QPE##2 

AXIS  C 1 , 10  » LOOP  > “AXIS  < 1 » 1 Or LOOP ) +Z  < 6 > #Z ( 8 ) 

AXIS  < 1 » 13  » LOOP ) =AXIS< 1 , 13 r LOOP > +QIE##2 
AXI S < 1 > 15 .LOOP )=AXIS<1, 15.L0QP ) +QAI#QAFI 

C 

C PASS  3 

C 

130  CONTINUE 

IF< IPASS  < NE . 3 ) GO  TO  140 
AX1S< 1 , 1 .LOOP )=AXIS< 1 . 1 .LOOP) +0A**2 
AXIS ( 1.3. LOOP) -AXIS< 1.3. L00P)+QE**2 
AXIS < 1 . 5. LOOP ) =AXIS< 1.5. LOOP ) +QA4QAF 
AXIS < 1 » 6.L00P ) -AXIS  < 1.6. LOOP ) +Z<  6 )#*2 
AXIS  < 1 . 7 .LOOP ) -AXIS  < 1 . 7 .LOOP ) +Z  < 8 > *#2 
AXIS  < 1 . 9 . LOOP ) -AXIS  <1.9. LOOP ) +0PE**2 
AXIS<  1. 11.  LOOP) -AXISU.il  .LOOP  >+Z<  6 >*Z<8> 

AXIS( 1 . 12. LOOP) -AX ISC  1 . 12,L00P)+QAI*#2 
AXIS <1,14. LOOP) -AXIS < 1 . 14 .LOOP ) +QIE**2  . 

AXIS  < 1 , 1 A > LOOP ) -AXIS  < 1 , 16  .LOOP ) +QAI*GAFI 

AXI S < 1 , 2.1 » LOOP ) - AM AXI  ( AXIS < 1 ,21  .LOOP ) , ABS < ZD ( 2 ) *57 . 3 > ) 

AXIS < 1, 22,  LOOP >-AMAXl< AXIS <1,22, LOOP > » ABS< Z< 2>*57 . 3) ) 

AXIS  < 1 , 23 , LOOP )-AXIS  < 1 , 23.L00P ) +0AF**2 

AXI S < 1 , 24 . LOOP ) -AXI 3 <1,24, LOOP ) +QAFI**2 

C 
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V 

I 


AX1 S<  2 » 1 . LOOP' ) - AXIS < 2.1.  LOOP ) +UP0D**2 
AXIS  < 2. 3 > LOOP ) =AXIS  < 2.3  .LOOF- ) +UED*#2 
AX IS  < 2 . 5 . LOOP ) =AX IS  < 2.  5 . LOOP  > +UPOD*UF'OFD 
AXI3<2.6.L00P>=AXIS<2.6.L00P>+Z<9>**2 
AXIS(2.7.L00P>=AXIS<2. 7. LOOP ) +Z< 12 >**2 
AXIS ( 2. 9 r LOOP ) -AXIS  < 2 . 9.L00P) +UPED442 
AXIS! 2 .11.  LOOP)  '-AXIS  ( 2.11  .LOOP )+Z  (9>#Z(12) 

AXIS ( 2 r 12  * LOOP ) “AXIS ( 2. 12  r LOOP ) +UP0ID##2 
AXIS!2»14.  LOOP  > “AXIS ( 2 » 1 4 » LOOP) +UIED**2 
AX I S ( 2 » 1 6 . LOOP ) “ AX I S ( 2 » 1 6 . LOOP  ) +UPO I D4UP0F I D 
AXIS <2f 23 » LOOP )=AXIS(2r23f LOOP )+UPOF»**2 
AXIS  < 2 . 24. LOOP > “AXIS ( 2 . 24 . LOOP >+UP0FID**2 
C 

C INTEGRATE 

C 

140  DO  150  1=1 .14 

Z< I)“Z0(I>+DELT2#(3.*ZD!I)-ZD0!I>) 
IF(I.EQ.6.0R.I.EQ.8>  Z ( I > = . 37*Z0( I >+ . 63*ZD( 1-1 > 
150  CONTINUE 
IT“IT+IDT 

IF(T.LT.TMAX)  GO  TO  90 
160  CONTINUE 
C 

C COMPUTE  MEANS 

C 

H- r MAX/DT 
DO  170  1=1.2 

DO  170  J=1 .24 

IF ( J.GE • 20. AND. J.LE .22)  GO  TO  170 
AXIS ( I . J. LOOP) “AXIS < I » J»L00P)/<N+1 > 

170  CONTINUE 

C 

IF (LOOP.NE • 1 ) GO  TO  200 
READ< IIN. 190)  NCOL 
190  FORMAT (II) 

200  L0GP=L00P+1 

IF  < LOOP • LE .NCOL > GO  TO  10 
C 

C PRINT  MEASURES 

XM  1 ) =XKQ 
XK (2  > “XK2P 
SCALE  < 1 > “XKQ 
SCALE! 2)=XK2P 
UNIT ( 1 ) "7.3 
UNIT<2)  0 
LUN<  1 ) = iiie.ASl 
LUN  < 2 ) “IMEAS2 
DO  228  1=1.8 
PEAK(1)«XKQ 
PEAK<2)=XK2P 
DO  228  J-1.24 
228  XMEAS  <I»J)“0.0 

c 

DO  270  IAX=1 .2 

IF ( SCALE ( I AX > <EQ .0.0)  SCALE < I AX ) = 1 .0 
IF" (PEAK ( IAX  > . EO. 0 • 0)  PFAK<  IAX)  = 1 .0 
DO  230  LOOP-1. NCOL 
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XMEAS (LOOPfI >=AXIS(IAXf 1 fLOOP )*UNIT( IAX >**2 

XMEA8 ( LOOP*  2 ) =AXXS  ( IAX  f2  f LOOP)/ AXIS ( I AX  f 1 r LOOP) /SCALE (IAX>*#2 
XMEAS(L00Pf4)=(XK(IAX)#AXIS(IAXf1  fLOOP  ) -AX  IS ( I AX  1 4 fl..  OOP  ) ) 

+ /AXIS(IAXf1fLOOP)/SCALE(IAX> 

XMEAS (LOOP » 3 ) =< AXIS(  I AX  > 3 . LOOP )-AXIS(  IAX  f2  fLOOP  > ) 

+ /AXISCIAXfI fL00P)/SCALE<IAX)**2 

XMEAS(L00Pf5)=(AXIS( IAXf4fL00P)~AXIS(IAXf5fL0QP) > 

+ /AXIS(IAXf1fL00P)/SCALE(IAX) 

XMEAS  ( L0QPf6)  - AXIS!  I AX  f 23  f LOOP ) #UNIT  ( I AX  > **2 
XMEAS (LOOP f 7 )=AXIS ( I AX  f6  fLOOP ) *UNIT ( l AX ) *42 
XMEAS ( LOOP  f 12 ) =AXIS ( I AX  f 7 fLOOP) #UNIT ( I AX ) **2 

XMEAS ( LOOP fB)=AXIS( IAX f8fL00P)/AXIS(IAXf6fL00P> /SCALE ( IAX >**2 
XMEAS  ( LOOP  f 10)  = ( XK  ( IAX  )*AXIS(  IAX  f6fL00P) -AXIS  (IAXf  10  fLOOP)  ) 

+ /AXIS(IAXf6fL00P)/SCALE(IAX) 

XMEAS  (LOOP'  f9  > - (AXIS  < I AXf  9 fLOOP  > - AXIS<  IAX  f8  fLOOP'  ) ) 

+ /AXIS ( IAX  f 6 fLOOP ) /SCALE ( I AX ) **2 

XMEAS (LOOP. 11 >=(AXIS( IAXf 10 fLOOP ) -AXIS ( IAXf 11 fLOOP) ) 

+ /AXIS(IAXf6fL00P)/SCALE(IAX> 

XMEAS (LOOP f 13 >=AXIS( IAX f12fL00P)#UNIT( IAX )**2 

XMEAS(  LOOP  f 14  >=AXIS(  IAXf  13fL.00P)/AXIS(  IAXf  12  fLOOP) /SCALE  ( IAX  >**2 
XMEAS (LOOP F 16>=(XK(IAX)*AXIS( IAXf 12fL00P)-AXIS( IAXf 15 fLOOP) ) 

+ /AXIS(IAXf12fL00P)/SCALE(IAX) 

XMEAS ( LOOP f 15 )=( AXIS! IAX f14fL00P)-AXIS(IAXf 13 fLOOP)) 

+ /AXIS(IAXf12fL00P)/SCALE(IAX>**2 

XMEAS ( LOOP  f17)  = (AXIS( IAXf 15  fLOOP ) -AXIS ( I AX  f 16  fLOOP) > 

+ /AXIS(IAXf12fL00P)/SCALE(IAX) 

XMEAS ( LOOP  f 18 ) =AXIS ( IAX  f 24  fLOOP ) *UNIT( IAX) **2 

XMEAS(L00Pf19>=0.0 

XMEAS(L00Pf20>=0.0 

XMEAS(L0OPf21)=O.O 

XMEAS ( LOOP f 22) =0.0 

XMEAS ( LOOP  f 23 ) =AX I S < I AX  f 21 f LOOP ) /PEAK (I AX ) 

XMEAS  < LOOP  f 24 ) = AX I S ( I AX  f 22  f LOOP ) /PEAK ( I AX ) 

IF< I AX • EQ . 2 ) XMEAS ( LOOP f 23 )=0«0 
IF ( IAX.EQ. 2)  XMEAS (LOOP f 24 > =0 • 0 
230  CONTINUE 
C 

DO  260  J=1f24 

WRITE ( LUN( I AX ) f 240  > ( XMEAS ( NLOOP  f J ) f NLOOP= 1 f NCOL ) 

240  FORMAT (8E) 

DO  245  NL00P=1fNC0L 
TEST=ABS(XMEAS(NLOOPfJ> ) 

IF(TEST. GE. 100.0. OR. (TEST. LE. 0.00005. AND. TEST. NE. 0.0))  GO  TO  252 
245  CONTINUE 

WRI TE  < I TERM  f 250 ) I L ABEL  < J ) f < XMEAS  < NLOOP  f J ) f NL00P=1 f NCOL ) 

250  FORMAT (1XfA5f1Xf8F9.5) 

GO  TO  260 

252  WRITE ( I TERM f 254)  ILABEL( J) f ( XMEAS (NLOOP fJ>  fNL00P=1 fNCOL) 

254  FORMAT ( IX  f A5f IX  f8E9 .2) 

260  CONTINUE 

WRITE ( I TERM f 265) 

265  FORMAT (///) 

270  CONTINUE 
C 

READ ( I IN  f 290 ) IANS 
290  FORMAT (13) 

IF( IANS.EQ.-l ) GO  TO  310 

C 

DO  300  1=1 f2 
DO  300  J*1f24 
DO  300  K=1f8 
300  AXIS(IfJfK)=0.0 
L00P=1 
GO  TO  10 
C 

310  STOP 
END 
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DIMENSION  Z<19)fZD<19>fZ0<l?)»ZD0(l?)»XK(2)rSCALE(2) 
DIMENSION  AXIS<2f24f8> fXMEASI 8f24 > f ILABEL < 24 > fUNIT < 2 > 
DIMENSION  LUN < 2 ) f XAUX < 8 > f Y < 19  > f YD  < 19  > 

DATA  IIN.ITERMfIMEAS1fIMEAS2  /20f5f21f22/ 

DATA  ( I LABEL (I>fI=1f24)  /5HPDLA  .5HPDL1  » 5HPDL2  f 
+ 5HPDL3  f 5HPDL4  f5HPDLS  i 3HPSENA 1 5HPSEN1 r 5HPSEN2  r 
+ 5HPSEN3 1 5HPSEN4 f 5HPSENS f 5HPNLA  fSHPNLI  f5HPNL2  r 5HPNL3  r 

+ 5MPNL4  t 5HPNLS  f5HGS1  f5HGS2  f5HGSS  fSHPACC  fdHPVEL  f 

+ SHPDIS  / 

LOOP=l 

C 

t:  GET  SCENARIO  FILENAME 

C 

10  READ( IIN.20)  NDEW 
20  FORMAT (12) 

C 

C CONSTANTS 

C 

READ ( I IN >30)  TMAXfDT 
,50  FORMAT  < 2F9 . 4 > 

0ELT2=DT/2. 

READ ( 1 1 N » 70 ) IORDER » UR  r C3Z f C4Z  f C5Z 
70  FORMAT (I1f4F9.4) 

C 

C INITIALIZE 

C 

DO  80  1=1 f 19 
Z<I>=0.0 
80  ZD< I >=0.0 

IT=0 
IDT=1 
C 

C DATA  COMPUTATION  LOOP 

C 

90  T=IT/10. 

DO  100  1=1 f 19 
ZO(I)=Z(I> 

100  ZDO< I )=ZD< I ) 

READ < NDEV f110>  PAfOAfRAf UPOD f VPOD f WPOD 
110  F0RMAT(A(F11.5f1X>) 

UP0D=WP0D/32.2+i. 

C 

C YAW 

C 

IF ( IORDER. EQ. 2)  GO  TO  120 
ZD(1)=0.5*RA-WR#Z<1) 

GO  TO  130 

120  ZD <1 >=0.5*RA-1.4*WR*Z(1)-UR**2*Z<2> 

ZD (2 *=Z<1> 

130  RAF=ZD< 1 ) 

ZD ( 3 ) “RA- . 2*Z ( 3 ) 

ZD(5)=RAF-.2*Z(5) 

RE=0.5*RA-RAF 
RPE=0.5*Z<4 )-Z<6) 

RAI=Z ( 4>-THRES< • 035 f Z < 4 > ) 

RAF I =Z ( 6 ) -THRES ( • 035  f Z ( 6 ) ) 

RIE=0.5*RAI-RAFI 


i 
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HEAtfF 

ZD  < 7 ) =0 . 15*WP0D- < . 05+C3Z+C5Z > *Z  < 7 >- < . 05*C5Z+C3Z*C5Z 
+ + . 05#C3Z+C4Z  > *Z  < 8 > - < . 05*C3Z*C5Z+C4Z*C5Z+ . 05*C4Z ) 

+ *Z<9>-.05*C4Z*C5Z*Z<10> 

ZD ( 8 ) =Z ( 7 ) 

ZD( 9 ) =Z  < 8) 

ZDIlO'ZfV) 

UP0FD=ZD<7>-1. 

ZD<11)-1.5*WP0D+.12*Z<12)-1.7*Z<11)-0.3*Z<13> 

ZD(12)=WP0D 

ZD<13)=Z<11> 

ZD<14)=1.5*<WP0FD+1.)+.12*Z<15>-1.7*Z<14>-0.3*Z<16) 

ZD ( 15 ) = < WP0FD+1 • > 

ZD< 1£)=Z< 14) 

ZD(17)=ZD(ll>+0.2*Z<lll-ra7> 

ZD< 18)=ZD( 14)+0.2*Z( 14)-Z(18) 

ZD<19)-0.92*WP0D-Z<19> 

WED-0. 15*UP0D-<UP0FD+1. ) 

UPED-0. 15#Z<  11  )-Z<  14) 

USED-0 .15*Z(17)~Z<18> 

W I ED-0. 15#<Z<17)+Z(19))-Z(18) 

UP0FSD»0.15*Z<19)-Z<18> 

AXIS< 1 r 1 f LOOP) -AX IS  < 1 f 1 fLOOP >+RA#*2 
AX I S < 1 f 3 f LOOP  > = AX IS  < 1 f 3 f LOOP ) +RE**2 
AXIS<1f5fL00P)«AXIS<1f5fL00P)+RA*RAF 
AXIS  < 1 f 6 f LOOP ) -AXIS  < l f 6 f LOOP  > +Z  < 4 )##2 
AXIS ( 1 f 7 f LOOP ) =AXIS< 1 f 7 f LOOP ) +Z ( 6 > #*2 
AXIS  < 1 f 9 fLOOP > -AXISZ1 f 9 f LOOP ) +RPE**2 
AXIS(1f11 fLOOP  > =AXIS( 1 f 1 1 f LOOP )+Z(4)*Z(6) 

AXIS<  1 f 12fLQ0F'>=AXIS<  1 f 12 fLOOP > +RAI**2 

AXIS ( 1 f 14  fLOOP ) “AXIS! 1 f 14fL00P ) +RIE**2 

AXIS< 1 f 16  fLOOP ) -AXIS< 1 f 16  fLOOP) +RAI4RAFI 

AXIS(1f21 fLOOP ) -AMAX1 < AXIS ( 1 f 21 fLOOP) f ABS( ZD  < 1 )$57 . 3 ) ) 

AXIS ( 1 f 22  fLOOP ) =AMAX1 ( AXIS( 1 f 22 fLOOP ) f ABS<Z ( 1 >$57 • 3 > ) 

AXIS<1f23 fLOOP )=AXIS( If  23 fLOOP >+RAF**2 

AXIS  < 1 f 24  f LOOP ) = AX I S < 1 f 24  f LOOP > +RAF I **2 

AXI S < 2 f 1 f LOOP ) -AX I S < 2 f 1 f LOOP  > +UP0D**2 

AXIS(2f3f  LOOP )-AXIS(2f3f  LOOP ) +UED442 

AXIS<  2f  5 f LOOP ) -AXIS  < 2 f5  f LOOP  > +UPOD*  < WPOFD+1 . ) 

AXIS ( 2 f6  fLOOP )-AXIS ( 2f  6fL00P) +Z< 1 1 >**2 
AXIS<  2 f 7 fLOOP ) -AXIS  < 2f  7fL00P) +Z  < 14 )#*2 
A X I S < 2 f 9 f L OOP ) - A X I S ( 2 f 9 f LOOP ) +UPED**2 
AXIS  ( 2 f 1 1 f LOOP  > “AXI  S<  2f  1 1 fLOOP)  +Z  ( 1 1 ) *Z  < 14  > 

AXIS(2f 12f LOOP)- AXI S<2f 12fL00P ) + <Z( 17)+Z< 19) )t$2 

AXIS(2f 14 fLOOP) -AXIS ( 2f 14fL00P)+USED**2 

AXIS( 2f 16 fLOOP) -AXIS <2f 16fL00P)+(Z< 17>+Z(19) > 4UP0FSD 

AXIS (2 f 17fL00P)-AXIS(2f 17 fLOOP >+Z< 18)4*2 

AXIS ( 2 f 18 fLOOP ) -AXIS<  2f 18 fLOOP ) +UIED442 

AXIS ( 2f 19fL00P ) -AXIS  < 2f 19 fLOOP ) + (Z<17)+Z<19) >*Z<18> 

AXIS(  2 f 20  fLOOP  ) - AM  AXI  < AXIS  < 2 f 20  fLOOP  ) f ADS  (ZD  (7X32  >2)) 

AXIS(2f21 fLOOP  > -AMAX1 (AXIS(2f21 fLOOP ) f APS  < Z ( 7 ) *32 . 2 > ) 

AX I S < 2 f 22 f LOOP ) -AMAX 1(AXIS(2f22f LOOP > f AB8 ( Z ( 8 ) *32 • 2 ) ) 

AXIS<2f23fLOOP  > -AXIS <2 f23 fLOOP > + < WPOFD+1 • >*#2 

AX IS(2f24fL00P> -AX I S ( 2 f 24  f LOOP ) +WP0FSD**2 
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C INTEGRATE 

C 

DO  140  1=1 f 19 

Z< I )=ZU( I )+DELT2*(3. *ZD( I >-ZDO< I ) > 

IF<I.E0.4.0R.I.EQ.6>  Z(  I ) = . 37*Z0( I > + . 63*ZD( 1-1 ) 

140  CONTINUE 
IT=IT+IHT 

IF(T.LT.TMAX)  GO  TO  90 
COMPUTE  MEANS 
N=TMAX/DT 

00  210  1=1 f 2 

DO  210  J=l»24 

IF < J . GE . 20 • AND < J . LE • 22 ) GO  TO  210 
AXIS! I. JfLOOP)=AXIS(If JfLOOP)/(N+l) 

10  CONTINUE 

IF<L00P.NE.1>  GO  TO  240 
READ! IINf 230)  NCOL 
30  FORMAT <I1) 

40  L00P=L00P+1 

IF(LOOP.LE.NCOL)  GO  TO  10 

PRINT  MEASURES 

XK( 1 )=0.5 
XK<2)=0« 15 
SCALE! 1 >=0.5 
SCALE (2  > =0< 15 
UNIT  < 1 )=57.3 
UNIT ( 2>=1 • 0 
LUN( 1 >=21 
LUN<2)=22 
DO  312  1 = 1 f 8 
DO  312  J=lr24 
312  XMEAS  (IfJ)=0.0 

C 

DO  360  IAX=1 >2 
DO  320  LOOP=lfNCOL 

XMEAS (LOOP » 1 ) "AXIS ( I AX  > 1 r LOOP > *UNIT ( IAX ) **2 

XMEAS  < LOOP  f 2 > -AXIS  < IAX , 3 1 LOOP > /AXIS  < I AX  > 1 r LOOP  > /SCALE  < I AX ) **2 
XMEAS (LOOP  f 4 > - ( XK ( I AX ) *AX IS ( I AX  * 1 r LOOP  > -AXIS ( I AX  1 5 1 LOOP  > ) 

4 /AXIS(IAXrl fLOOP ) /SCALE ( IAX ) 

XMEAS ( LOOP  * 6 ) -AXIS ( IAX , 23 . LOOP ) *UNIT ( I AX ) **2 
XMEAS ( LOOP  > 7 ) -AX IS(IAXr6r  LOOP  > *UNI T ( I AX ) **2 

XMEAS (LOOP  1 8 ) -AXIS ( I AX  > 9 > LOOP ) /AXIS ( IAX  r 6 fLOOP  > /SCALE ( I AX ) **2 
XMEAS ( LOOP f 10>-(XK( IAX ) *AXIS( IAX? 6f LOOP ) -AXIS ( IAXr 1 1 r LOOP) ) 

+ /AXIS( I AXr6f LOOP) /SCALE! IAX) 

XMEAS ( LOOP  r 1 2 ) -AXIS ( I AX  r 7f LOOP ) *UNIT ( I AX  > **2 
XMEAS (LOOP  r 13) -AXIS ( IAX  1 12f  LOOP ) *UNIT( I AX ) **2 

XMEAS(LOOPf 14 > -AXIS (IAXf 14 fLOOP >/AXIS( I AXfl2f LOOP) /SCALE (IAX) **2 
XMEAS (LOOPf 16 )-( XK< I AX  >*AXIS( IAXf 12f LOOP) -AXIS! IAXf 16f LOOP) ) 

+ /AXIS ( IAXf 12f LOOP )/SC ALE ( I AX ) 

XMEAS  ( LOOP  f 18  > =AX  IS  ( I AX  f 24  fLOOP  > *UNI  T ( I AX ) **2 

IF ( I AX .EG. 1 ) 00  TO  315  , 

XMEAS ( LOOP  f 19  > -AXIS ( I AX  f 18f  LOOP ) /AX IS ( IAX  f 1 2 f LOOP ) /SCALE ( I AX  > **2 
XMEAS( LOOPf 20 >=(XK( IAX )*AXIS( IAXf 12f LOOP >-AXIS( IAXf 19 fLOOP)) 
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+ /AXIS< IAXf 12 » LOOP) /SCALE ( I AX) 

XMEAS < LOOP f 21 ) =AX IS ! I AX  f 1 7 f LOOP  > *UNIT ( I AX ) **2 
XMEAS! LOOP f 22 )=AXIS! IAXf  20  r LOOP  > /XK  < I AX ) 

315  XMEAS  < LOOPf  23 )=AXIS! I AX  >21 r LOOP) /XK  < I AX ) 

XMEAS < LOOP*  24  > “AXIS! IAXf 22 » LOOP ) /XK! IAX  > 

320  CONTINUE 
C 

DO  350  J-l i 24 

WRI TE ! LUN ! I AX  > f 325  > < XMEAS  < NLOOP  f J > f NI  OOP=  1 f NCOL ) 

325  FORMAT !8E> 

DO  330  NL00P=1fNC0L 
TEST= ADS ! XMEAS  < NLOOP f J > > 

IF ! TEST. GE . 100.0. OR. < TEST* LE. 0 . 00005. AND. TEST . NE .0 .0) ) GO  TO  340 
330  CONTINUE 

WRI TE ! I TERM  f 335 ) t LABEL ( J ) t ( XMEAS ! NLOOP  fJ)f NLOOP= 1 f NCOL > 

335  FORMAT ( IXfASf 1Xf8F9.5) 

GO  TO  350 

340  WRITE! I TERM f 345)  ILABEL < J) f (XMEAS (NLOOP fJ ) fNLOOP-1 fNCOL  ) 

345  FORMAT (1XfA5f1Xf8E9.2) 

350  CONTINUE 

WRITE! I TERM f 355) 

355  FORMAT!///) 

360  CONTINUE 
C 

READ! I IN f 380)  IANS 
380  FORMAT (13) 

IF! IANS.EQ. -1 ) GO  TO  400 
C 

DO  390  1=1 f2 
DO  390  J=1 f24 
DO  390  K=1f8 
390  AXIS!IfJfK)=0.0 
L00P=1 
GO  TO  10 
C 

400  RETURN 
END 
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C****THRES-START 

FUNCTION  THRESITHfX) 

C 

C IMPLEMENTS  A THRESHOLD.  OF  THE  FORM 

C y = TH*SIN<X/TH)f  FOR  -PI#TH/2<X<PI*TH/2 

C AND 

C Y = -TH  FOR  Xf~PI*TH/2 

C T=TH  FOR  X>PI*TH/2 

C 

DPR  = 57.29578 

PI  = 3.14159 

IFIX.LT .0.0)  THRES=-TH 

IF(X.GT.O.O)  THRES=TH 

IF ( X . GE . -PI* TH/2. . AND . X .LE . PI*TH/2. ) 

+ THRES  = TH*SIN(X/TH) 

RETURN 

END 

C*«**THRES-END 

C*$**THRESD-START 

FUNCTION  THRESD(THfX) 

C 

C IMPLEMENTS  A THRESHOLD  OF  THE  FORM 

C Y=C0S<PI*X/2/TH)f  FOR  -TH<X<TH 

C AND 

C Y-0  f FOR  X<-TH  OR  X>TH 

C 

DPR*57. 29578 
Pl=3. 14159 
THRESD=0 • 0 

I F < X . GE . -TH . AND . X . LE . TH ) THRESD=COS ( X/TH ) 

RETURN 

END 


C****THRESD-END 


SECTION  III 


INPUT  AND  OUTPUT 


Most  of  the  input  and  output  involved  in  this  subroutine  package 
is  stored  on  files  — the  only  exceptions  to  this  are  the  hard  copy 
plots  which  can  be  produced,  and  the  listing  of  measures  via  the  line 
printer.  In  this  section,  each  subroutine  which  requires  input  or  pro- 
duces output  will  be  discussed  in  terms  of  detailed  descriptions  of  that 
input  or  output.  At  the  end  of  this  section  are  included  sample  listings 
of  input  and  output. 

MAINLINE 

For  the  scenario  generation  mainline,  all  of  the  input  data  is  con- 
tained on  the  logical  device  IIN.  The  first  input  is  the  list  of  possible 
plotting  variable  names.  There  must  be  150  names  input:  50  state  names, 

50  state  derivative  names  and  50  auxiliary  variable  names,  in  that  order. 
Each  name  is  5 characters  long.  The  file  must  contain  25  lines,  6 names 
to  a line,  separated  by  1 comma  or  1 blank.  All  names  with  fewer  than  5 
characters  must  be  filled  out  with  blanks.  Each  state  variable  name  must 
correspond  to  its  appropriate  state  — that  is,  the  first  state  variable 
name  [NVLA  (1)]  must  correspond  to  the  first  state  [Z ( 1 ) ] • The  same 
holds  true  for  the  state  derivatives  and  auxiliary  variables  [NVLA(51 ) 
corresponds  to  ZD(l)  and  NVLA(ioi)  corresponds  to  XAUX(l)].  If  there 
are  fewer  than  50  states,  state  derivatives  or  auxiliary  variables,  5 
blanks  should  comprise  the  remaining  variable  names.  The  next  input  is 
the  number  of  variables  to  be  plotted.  This  is  an  integer  greater  than  0 
and  less  than  21 . Next  are  the  plotting  variable  names  and  their  corres- 
ponding scale  factors.  Each  variable  name  and  its  scale  factor  must 
occupy  a separate  line.  The  variable  name  must  correspond  to  a name  in 
the  list  of  possible  variable  names.  The  variable  name  must  be  5 charac- 
ters long,  the  scale  factor  is  in  units  per  inch  and  must  be  a real  number, 
and  the  two  must  be  separated  by  a comma  or  a blank.  The  final  input  used 
by  this  routine  is  the  maximum  time,  in  seconds  and  the  time  scale  factor, 
in  seconds  per  inch.  Each  is  a real  number. 
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The  output  from  this  subroutine  is  a set  of  hard  copy  plots.  Up 
to  20  plots  are  permitted,  five  per  page*  These  plots  are  output  to 
device  IPLOT* 

SXHGKN 

There  are  several  inputs  to  this  subroutine  — they  should  follow 
directly  after  the  data  for  the  mainline  discussed  above.  The  first  is 
the  scenario  number,  NSCEN.  It  is  an  integer  which  uniquely  identifies 
which  of  the  eight  existing  scenarios  is  desired.  Hie  next  line  con- 
tains an  integer  whose  value  is  the  number  of  discrete  data  points  needed 
to  describe  this  scenario  (NDISC).  A list  of  the  scenarios  and  their 
corresponding  NSCEN  and  NDISC  values  is  given  below: 


Scenario 

NSCEN 

NDISC 

HQPTl 

1 

5 

HQDT2 

2 

9 

HQDT4 

3 

- 

AG 

4 

4o 

AA1 

5 

42 

AA2 

6 

16 

AA3 

7 

8 

AA4 

8 

9 

The  final  input  for  scenario  generation  is  a list  of  the  discrete  values 
from  which  the  time  history  of  the  scenario  will  be  computed.  (Note, 
however,  that  the  HQDT4  scenairio  does  not  require  the  input  of  NDISC  or 
the  list  of  discrete  values.  This  is  because  the  HQDT4  sceneirio  is 
described  by  algebraic  equations  which  are  functions  of  time,  not  by  dis- 
crete values).  The  discrete  values  include  TMAX,  PA,  QA  and  VTDD-  The 
TMAX  value  is  the  maximum  value  for  which  the  PD,  QD  and  VTDD  values 
hold  — when  one  or  more  of  these  values  chatnges,  a new  TMAX  is  needed. 


A more  complete  description  of  the  process  used  to  arrive  at  the  dis- 
crete values  can  he  found  in  the  Scenario  Generation  chapter  of  the 
Engineering  Flight  Simulation  Validation  final  report. 

The  output  from  the  scenario  generation  is  a file  or  device  IOUT 
which  contains  the  inputs  to  the  drive  logic  simulation  (PA,  QA,  RA, 
UPOD,  VPQD,  WPQD) . Each  of  these  variables  is  computed  at  0.1  second 
intervals  for  the  length  of  the  scenario  (TMAX) . 

ROLSWA,  PITSUR,  YAWHEV 


The  inputs  to  these  routines  are  all  very  similar,  and  will  be  dis- 
cussed together.  Tht.  a are  actually  two  sets  of  inputs  — the  first  is 
the  scenario,  which  resides  on  logical  unit  NDEV.  This  file  was  pro- 
duced by  the  scenario  generation  package.  The  other  inputs  reside  on 
logical  unit  IIN.  These  inputs  are  described  below: 


NDEV 

TMAX,DT 

Constants 

For  ROLSWA 

For  PITSUR 
For  YAWHEV 

NCOL 

IANS 


logical  unit  for  scenario  generation  data 

maximum  length  of  scenario,  and  integration 
interval 

(for  definition  of  constants,  see  Final 
Report) 

KP,  WP,  W1R,  W2R,  W5R,  W^R,  K2R,  C2Y 
Drive  logic  type  (l,  2 or  5),  BP,  BK,  C2YP 
K2P,  WQ,  W3P,  W^P 
IORDER  (l  or  2),  WR,  C3Z,  CUZ,  C5Z 

number  of  columns-  If  more  than  1 column 
of  data  is  desired,  all  of  the  preceding 
information  must  be  supplied  for  each  column. 

If  more  data  is  to  be  computed,  a 0 is  input 
and  all  preceding  information  must  be  included 
for  the  next  page  of  data-  Otherwise,  a — i 
input  ends  the  program- 


Outputs  from  each  of  the  three  routines  are  files  containing  the 
angular  velocity  and  specific  force  axes  measures  (on  devices  IMEAS1 
and  IMEAS2,  respectively),  as  well  as  a hard  copy  formatted  listing  of 
those  measures. 


SECTION  17 


A TYPICAL  APPLICATION 


This  example  presents  the  input  files,  output  files  and  hard  copy 
associated  with  generating  the  AA4  scenario  and,  subsequently,  the  per- 
formance measures  for  two  different  sets  of  drive  logic  parameters  when 
the  drive  logic  forcing  function  is  the  AA4  scenario.  Only  the  roll  and 
sway  measures  are  presented.  Each  section  of  each  file  is  labeled  to 
indicate  what  variable  or  set  of  variables  is  listed. 

A listing  of  the  discrete  values  used  to  generate  each  of  the  seven 
remaining  scenarios  follows  this  example. 


LIST  OF  INPUTS  AND  OUTPUTS  FOR 
SAMPLE  RUN  (Scenario  8) 


TINFUT  — Contains  all  inputs  to  scenario  generation 
software . 

Scenario  generation  plots  for  Scenario  8 

TOUT  — Contains  time  histories  for  PA,  QA,  RA,  UPOD, 

VPOD,  WPOD  variables.  This  file  is  output  by 
the  scenario  generation  software,  and  is  read 
as  input  to  the  drive  logic/sensory  apparatus 
simulation  software. 

DINPUT  - Contains  all  inputs  to  drive  logic/sensory  appa- 
ratus simulation  software  except  those  contained 
in  TOTJT 

Table  of  measures  generated  by  drive  logic/sensory  apparatus 

simulation,  for  roll- sway  axes. 

ROLL  — File  containing  roll  axis  measures 

SWAY  — File  containing  sway  axis  measures 

I 
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1?  - MV 
PD  ,50. 

PA  .50. 
PHID  .50. 

PHI  .100. 
VPOD  ,.5 
QD  .50. 

QA  .50. 

THD  .50. 

TH  .50. 
UPOD  ,2. 

VT  .2000. 
RA  .50. 
PSID  .50. 

PSI  .360. 
UPOD  .10. 

VTD  ,100. 
VTDD  ,50. 
ALFA  ,50. 
ALFAD.50. 

15.0  ‘ 

8 - NSCEN 

9 - NDISC 

2.0 
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RD 

ALFAD 


- KVLA,  SF 


5.0  - TMAX,  TIMSF 


-.192 

.314 

0.0 

-.960 

1.030 

.175 

-.175 

0.0 

0.0 


0.0 

0.0 

0.0 

.384 

0.0 

-.349 

-.349 

.349 

0.0 


- TMAX,  PD,  Sp,  VTDD 


I 


> 
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0.00000  0.00000 

-0.02880  0.00000 

-0.04800  0.00000 

-0 . 06720  0.00000 

-0.08640  0.00000 

-0.10560  0.00000 

-0.12480  0.00000 

-0.14400  0.00000 

-0.16320  0.00000 

-0.18240  0.00000 

-0.20160  0.00000 

-0.22080  0.00000 

-0.24000  0.00000 

-0.25920  0.00000 

-0.27840  0.00000 

-0.29760  0.00000 

-0.31680  0.00000 

-0.33600  0.00000 

-0.35520  0.00000 

-0.37440  0.00000 

-0.39360  0.00000 

-0.33690  0.00000 

-0.30550  0.00000 

-0.27410  0.00000 

-0.24270  0.00000 

-0.21130  0.00000 

-0.17990  0.00000 

-0.14850  0.00000 

-0.11710  0.00000 

-0.08570  0.00000 

-0.05430  0.00000 

-0.07000  0.00000 

-0.07000  0.00000 

-0.07000  0.00000 

-0.07000  0.00000 

-0.07000  0.00000 

-0.07000  0.00000 

-0.07000  0.00000 

-0.07000  0.00000 

-0.07000  0.00000 

-0.07000  0.00000 

-0.21400  0.05760 

-0.31000  0.09600 

-0.40600  0.13440 

-0.50200  0.17280 

-0.59800  0.21120 

-0.69400  0.24960 

-0.79000  0.28800 

-0.88600  0.32640 

-0.98200  0.36480 

-1.07800  0.40320 

-0.87550  0.38400 

-0.77250  0.38400 

-0.66950  0.38400 

-0.56650  0.38400 

-0.46350  0.38400 

-0.36050  0.38400 

-0.25750  0.38400 
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0.00000  2.25400 
-0.00202  2.25375 
-0.00353  2.25327 
-0.00511  2.25251 
-0.00676  2.25142 
-0.00849  2.24993 
-0.01030  2.24796 
-0.01218  2.24538 
-0.01414  2.24206 
-0.01618  2.23782 
-0.01829  2.23243 
-0.02047  2.22564 
-0.02273  2.21714 
-0.02505  2.20658 
-0.02744  2.19356 
-0.02989  2.17764 
-0.03240  2.15830 
-0.03497  2.13501 
-0.03758  2.10716 
-0.04023  2.07412 
-0.04291  2.03520 
-0.04040  2.01380 
-0.03926  1.97763 
-0.03806  1.93798 
-0.03674  1.89389 
-0.03532  1.84614 
-0.03381  1.79531 
-0.03220  1.74205 
-0.03050  1.68702 
-0.02872  1.63087 
-0.02686  1.57426 
-0.02799  1.51063 
-0.02823  1.45062 
-0.02843  1.38973 
-0.02864  1.32878 
-0.02885  1.26755 
-0.02905  1.20598 
-0.02926  1.14399 
-0.02947  1.08155 
-0.02967  1.01858 
-0.02987  -0.23534 
-0.03914  -6.38861 
-0.04870  -9.86613 
-0.06069  -13.11797 
-0.07618  -15.94567 
-0.09568  -18.26002 
-0.11965  -19.96976 
-0.14845  -20.99700 
-0.18234  -21.28271 
-0.22153  -20.79085 
-0.26615  -18.32043 
-0.24625  -8.24904 
-0.23570  -2.73634 
-0.22114  2.92372 
-0.20204  8.14367 
-0.17915  12.93039 
-0.15305  17.27410 
-0.12423  21.18317 


-0.59520  -32.20000 

-0.59520  -32.19640 

-0.59520  -32.18986 

-0.59520  -32.17989 

-0.59520  -32.16627 

-0.59520  -32.14867 

-0.59520  -32.12666 

-0.59520  -32.09963 

-0.59520  -32.06682 

-0.59520  -32.02725 

-0.59520  -31.97978 

-0.59520  -31.92304 

-0.59520  -31.85542 

-0.59520  -31.77509 

-0.59520  -31.67998 

-0.59520  -31.56779 

-0.59520  -31.43595 

-0.59520  -31.28167 

-0.59520  -31.10189 

-0.59520  -30.89335 

0.97340  -30.65256 

0.97340  -30.64447 

0.97340  -30.50544 

0.97340  -30.34492 

0.97340  -30.15900 

0.97340  -29.95352 

0.97340  -29.73301 

0.97340  -29.50218 

0.97340  -29.26575 

0.97340  -29.02830 

0.00000  -28.79429 

0.00000  -28.52980 

0.00000  -28.29978 

0.00000  -28.07338 

0.00000  -27.85677 

0.00000  -27.64781 

0.00000  -27.44570 

0.00000  -27.24968 

0.00000  -27.05905 

0.00000  -26.87318 

-2.97600  -33.47679 

-3.18682  -33.42867 

-3.48737  -40.20672 

-3.91490  -48.11546 

-4.46799  -58.06284 

-5.14503  -69.71411 

-5.94426  -82.81303 

-6.86375  -97.10711 

-7.90150  -112.35280 

-9.05547  -128.31628 

-4.15460  -137.98941 

-2.45438  -158.12007 

-1.76805  -171.07873 

-1.08650  -183.33646 

-0.41035  -193.93938 

0.26128  -203.12115 

0.92909  -211.02662 

1.59367  -217.79003 
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-0.15450 

0.38400 

-0.09310 

24.47738 

2.25550 

-223.53195 

-0.05150 

0.38400 

-0.05999 

27.78449 

2.91497 

-228.34071 

0.05150 

0.38400 

-0.02514 

31.41952 

0.92190 

-226.20664 

0.02425 

0.33145 

-0.03398 

34.70709 

0.73124 

-229.72482 

0.04375 

0.29475 

-0.02810 

33.73442 

0.79775 

-227.20657 

0.04125 

0.24185 

-0.02240 

31.80390 

0.84408 

-223.25377 

0.07875 

0.22495 

-0.01705 

28.57249 

0.87030 

-217.24677 

0.09425 

0.19205 

-0.01218 

24.27924 

0.87438 

-209.52364 

0.11375 

0.15715 

-0.00789 

19.13745 

0.86221 

-200.23588 

0.13125 

0.12225 

-0.00429 

13.35412 

0.82745 

-189.59441 

0.14875 

0.08735 

-0.00147 

7.13494 

0.77249 

-177.77011 

0.14425 

0.05245 

0.00050 

0.44170 

0.69650 

-164.91599 

0.18375 

0.01755 

0.00155 

-5.88924 

-0.48560 

-151.16913 

0.14875 

-0.01735 

-0.01089 

-4.85384 

-0.58869 

-134.95767 

0.13125 

-0.05225 

-0.01817 

-4.04045 

-0.66662 

-119.19783 

0.11375 

-0.08715 

-0.02482 

-4.90771 

-0.72437 

-103.05099 

0.09425 

-0.12205 

-0.03073 

-7.33338 

-0.76170 

-86.60418 

0.07875 

-0.15495 

-0.03583 

-7.21423 

-0.77833 

-69.87564 

0.04125 

-0.19185 

-0.04005 

-4.45414 

-0.77396 

-52.87619 

0.04375 

-0.22475 

-0.04333 

-4.95428 

-0.74825 

-35.40723 

0.02425 

-0.24145 

-0.04541 

-2.41187 

-0.70086 

-18.06194 

0.00875 

-0.29455 

-0.04483 

0.47048 

-0.63140 

-0.22595 

-0.00875 

-0.33145 

-0.04494 

2.83091 

0.00300 

5.58817 

-0.00000 

-0.24145 

-0.04490 

1.94415 

-0.03804 

24.11995 

-0.00000 

-0.22475 

-0.04484 

3.14809 

-0.03292 

34.02487 

-0.00000 

-0.19185 

-0.04477 

4.44108 

-0.02781 

42.51945 

-0.00000 

-0.15495 

-0.04470 

5.44271 

-0.02272 

48.36358 

-0.00000 

-0.12205 

-0.04443 

4.04140 

-0.01764 

51.81780 

-0.00000 

-0.08715 

-0.04457 

4.24855 

-0.01258 

53.05794 

-0.00000 

-0.05225 

-0.04450 

4.00990 

-0.00753 

52.24964 

-0.00000 

-0.01735 

-0.04443 

5.39385 

-0.00250 

49.54589 

-0.00000 

0.01755 

-0.04434 

4.47974 

0.00252 

45.08815 

-0.00000 

0.05245 

-0.04430 

4.45131 

0.00753 

45.17400 

-0.00000 

0.03500 

-0.04423 

3.37520 

0.00502 

37.57075 

-0.00000 

0.03500 

-0.04414 

2.77290 

0.00501 

33.27074 

-0.00000 

0.03500 

-0.04410 

2.21299 

0.00500 

28.65718 

-0.00000 

0.03500 

-0.04403 

1.75444 

0.00499 

24.43748 

-0.00000 

0.03500 

-0.04597 

1.37474 

0.00499 

20.54844 

-0.00000 

0.03500 

-0.04590 

1.04443 

0.00498 

16.96583 

-0.00000 

0.03500 

-0.04584 

0.80529 

0.00497 

13.66563 

-0.00000 

0.03500 

-0.04577 

0.58897 

0.00497 

10.62574 

-0.00000 

0.03500 

-0.04571 

0.40489 

0.00496 

7.82576 

-0.00000 

0.03500 

-0.04544 

0.25195 

0.00495 

5.24686 

-0.00000 

0.03500 

-0.04558 

0.11827 

0.00495 

2.87165 

-0.00000 

0.03500 

-0.04551 

0.00105 

0.00494 

0.68410 

-0.00000 

0.03500 

-0.04545 

-0.10349 

0.00493 

-1.33059 

-0.00000 

0.03500 

-0.04538 

-0.19914 

0.00492 

-3.18608 

-0.00000 

0.03500 

-0.04532 

-0.28800 

0.00492 

-4.89497 

-0.00000 

0.03500 

-0.04524 

-0.37232 

0.00491 

-6.46886 

-0.00000 

0.03500 

-0.04519 

-0.45382 

0.00490 

-7.91847 

-0.00000 

0.03500 

-0.04513 

-0.53387 

0.00490 

-9.25368 

-0.00000 

0.03500 

-0.04507 

-0.41352 

0.00489 

-10.48359 

-0.00000 

0.03500 

-0.04500 

-0.49341 

0.00488 

-11.61660 

-0.00000 

0.03500 

-0.04494 

-0.77477 

0.00488 

-12.66045 

-0.00000 

0.03500 

-0.04488 

-0.85749 

0.00487 

-13.62226 

-0.00000 

0.03500 

-0.04482 

-0.94212 

0.00486 

-14.50861 

-0.00000 

0.03500 

-0.04475 

-1.02889 

0.00486 

-15.32555 

-0.00000 

0.03500 

-0.04449 

-1.11797 

0.00485 

-16.07865 
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-0.00000 

0.03500 

-0.04463 

-1.20945 

0.00484 

-16.77304 

-0.00000 

0.03500 

-0.04457 

-1.30337 

0.00484 

-17.41347 

-0.00000 

0.03500 

-0.04451 

-1.39972 

0.00483 

-18.00428 

-0.00000 

0.03500 

-0.04444 

-1.49847 

0.00482 

-18. 54949 

-0.00000 

0.03500 

-0.04438 

-1.59955 

0.00482 

-19.05278 

-0.00000 

0.03500 

-0.04432 

-1.70288 

0.00481 

-19.51755 

-0.00000 

0.03500 

-0.04426 

-1.80837 

0.00480 

-19.94694 

-0.00000 

0.03500 

-0.04420 

-1.91591 

0.00480 

-20.34380 

-0.00000 

0.03500 

-0.04414 

-2.02540 

0.00479 

-20 .71081 

-0.00000 

0.03500 

-0.04408 

-2.13673 

0.00478 

-21.05038 

-0.00000 

0.03500 

-0.04402 

-2.24978 

0.00478 

-21.36475 

-0.00000 

0.03500 

-0.04396 

-2.36444 

0.00477 

-21.65600 

-0.00000 

0.03500 

-0.04390 

-2.48061 

0.00476 

-21.92601 

-0.00000 

0.03500 

-0.04384 

-2.59816 

0.00476 

-22.17653 

-0.00000 

0.03500 

-0.04378 

-2.71701 

0.00475 

-22.40915 

-0.00000 

0.03500 

-0.04372 

-2.83704 

0.00474 

-22.62535 

-0.00000 

0.03500 

-0.04366 

-2.95818 

0.00474 

-22.82648 

-0.00000 

0.03500 

-0.04360 

-3.08032 

0.00473 

-23.01378 

-0.00000 

0.03500 

-0.04354 

-3.20338 

0.00472 

-23.18840 

-0.00000 

0.03500 

-0.04348 

-3.32728 

0.00472 

-23.35138 

-0.00000 

0.03500 

-0.04343 

-3.45195 

0.00471 

-23.50370 

-0.00000 

0.03500 

-0.04337 

-3.57731 

0.00471 

-23.64622 

-0.00000 

0.03500 

-0.04331 

-3.70331 

0.00470 

-23.77978 

-0.00000 

0.03500 

-0.04325 

-3.82987 

0.00469 

-23.90511 

-0.00000 

0.03500 

-0.04319 

-3.95694 

0.00469 

-24.02290 

-0.00000 

0.03500 

-0.04313 

-4.08447 

0.00468 

-24.13377 

-0.00000 

0.03500 

-0.04308 

-4.21240 

0.00467 

-24.23832 

-0.00000 

0.03500 

-0.04302 

-4.34070 

0.00467 

-24.33705 

-0.00000 

0.03500 

-0.04296 

-4.46932 

0.00466 

-24.43047 

-0.00000 

0.03500 

-0.04290 

-4.59821 

0.00466 

-24.51902 

-0.00000 

0.03500 

-0.04285 

-4.72735 

0.00465 

-24.60311 

-0.00000 

0.03500 

-0.04279 

-4.85670 

0.00464 

-24.68310 

-0.00000 

0.03500 

-0.04273 

-4.98622 

0.00464 

-24.75936 

-0.00000 

0.03500 

-0.04268 

-5.11590 

0.00463 

-24.83218 

-0.00000 

0.03500 

-0.04262 

-5.24569 

0.00462 

-24.90187 

-0.00000 

0.03500 

-0.04256 

-5.37559 

0.00462 

-24.96869 
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SUBLIMINAL  LINEAR 
('OLA  J.20E+03  U.23E+03  > 

P0L1  1.03192  1.03192  j 

HOL2  0.63262  0.63262 

('DL3  0,70101  3.70101 

(’UL4  0. 3340  3 0.33483 

POLS  0.22039  0.22039 

('SENA  0.16E+0.J  0.16E+03 

PSENl  1.02612  1.02612 

(’SEN2  0.73413  0.73413 

('SEN3  0. 726^7  0.72697 

(’SEN4  0.4)207  0.43287 

(’SENS  0.23436  0.20436 

(’NLA  0.13E+33  0.13E*03  Jp  measures 

(’NL1  0.96801  0.93001 

(’NL2  0.04672  0.04672 

l>flL3  0.98307  0.90337 

PNL4  3.32334  0.02334 

(•NLS  0.00040  0.33348 

081  0.00103  0.03300 

QS2  0.00300  0.30003 

GSS  0.00300  0.03303 

('ACC  0.00003  0.00030 

(’VEL  33.61460  30.61460 

(’DIS  29.7?'’43  2a  .7Q,’46  , 


(’OLA  3.00266  0.00266 

POL1  0.68963  3.60925 

H0L2  0.33173  3.00381 

POL3  0.34229  0.3039’ 

PDL4  -3.32917  -3.33374 
COLS  0.33316  0.30310 

(’SENA  3.33133  0.30130 

(’SEMI  0.55926  0.33633 

(’5EN2  1. Oll^O  3.01635 

(>SEN3  0.36271  0.33645 

PSFN4  -0.33837  0.31462 

(’SENS  0.78E-04  3.37E-34  Jy  measures 

(’NLA  0.48E— 34  0.48E-04  T 

PNL1  0.965’5  3.92194  f 

r NL2  -0.002'’2  0.33253 

PNL3  0.98233  0.96313 

PNL4  -0.03165  3.00132 

I’Nt.S  0.13E-03  0.27E-08  ' 

GSl  3.30030  3.33333 

GS2  3.33303  0.33033 

GS'j  0.03301  3,30113 

PAC'"  27.9791?  36.29805 
PVEL  16.43408  24.66418 
POIS  21.68047  62.36570 


POLA  3.00251  0.00281  \ 

POL l 0.38290  0.00414  I 

(’01.7  0.59564  0.00901  > J measures 

POL3  0.24131  0.02167  I ' 

P0L4  -0.06870  0.30157  I 
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ROLL 


MON  31-JUL-78  13103 


203.65 

203.65 

1.03 

1.03 

0.63 

0.63 

0.78 

0.78 

0.33 

0.33 

0.22 

0.22 

157.18 

157.18 

1.03 

1 >03 

0.75 

0.75 

0.73 

0.73 

0.40 

0.40 

0.20 

0.20 

126.16 

126.16 

0.97 

0.97 

0.05 

0.05 

0.98 

0.98 

0.02 

0.02 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

30.61 

30.61 

29.80 

29.80 
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SWAY 


MON  31-JUL-78  13J04 


0.00 

0.00 

0.69 

0.61 

0.00 

0.00 

0.84 

0.80 

-0 . 03 

-0.03 

0.00 

0.00 

0.00 

0.00 

0.56 

0.01 

1.08 

0.02 

0.36 

0,04 

-0.01 

0.01 

0.00 

0.00 

0.00 

0.00 

0.97 

0.92 

-0.00 

0.00 

0.98 

0.96 

-0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

27.98 

36.30 

16.43 

24,66 

21.68 

52.07 

0.00 

0.00 

0.38 

0.00 

0.60 

0.01 

0.24 

0.02 

-0.07 

0.00 

a 


DISCRETE  DATA  FOR  REMAINING  SCENARIOS 


Scenario  1 , HQDTl : 5 discrete  data  points 


2 ♦ 0000 

0.0000 

0,0000 

0.0000 

3*0000 

0.4400 

0.0000 

0.0000 

5.0000 

0.0000 

0.0000 

0.0000 

6.0000 

-0.4400 

0.0000 

0.0000 

11.0000 

0.0000 

0.0000 

0.0000 

Scenario  2,  HQDT2:  9 discrete  data  points 

2.0000 

0.0000 

0.0000 

0.0000 

3.0000 

0.4400 

0.0000 

0.0000 

5.0000 

0.0000 

0.0000 

0.0000 

6.0000 

-0.4400 

0.0000 

0.0000 

8.0000 

0.0000 

0.0000 

0.0000 

9 . 0000 

-0.6600 

0,0000 

0.0000 

12.0000 

0.0000 

0.0000 

0.0000 

13.0000 

0.6600 

0.0000 

0.0000 

17.0000 

0.0000 

0.0000 

0.0000 

Scenario  4,  A/G: 

4o  discrete 

data  points 

3.0000 

0.0000 

0.0000 

0.0000 

5.0000 

0.0000 

0.0000 

2.5000 

6,0000 

-0.5200 

0.0350 

0.0000 

7.0000 

0.0000 

0,0350 

0.0000 

8.0000 

0.5200 

0.0350 

0.0000 

22.4000 

0.0000 

0.0000 

0.0000 

23.4000 

0.5200 

-0.0350 

0.0000 

24,4000 

0.0000 

-0.0350 

0.0000 

25.4000 

-0.5200 

-0.0350 

0.0000 

45.0000 

0.0000 

0.0000 

0.0000 

46.0000 

-0.8700 

0.0750 

4.0000 

47 . 0000 

0.0000 

0.0750 

4.0000 

48.0000 

0.8700 

0.0750 

4.0000 

52.0000 

0.0000 

0.0000 

0.0000 

53.0000 

0.8700 

-0.0750 

0.0000 

54.0000 

0.0000 

-0.0750 

0.0000 

55.0000 

-0.8700  - 

-0.0750 

0.0000 

65.5000 

0.0000  - 

0.0000 

0.0000 

66.5000 

0.0000 

0.1640 

-4,0000 

67.0000 

0.0000 

0.0000 

-4.0000 

68.0000 

-0.7000 

0.0070 

-4.0000 

68.5000 

0.0000 

0.0070 

-4.0000 

69.0000 

0.0000 

0.0070 

-2.5000 

70.0000 

0.7000 

0.0070 

-2.5000 

70.5000 

0.0000 

0.0070 

-2.5000 

DISCRETE  DATA  FOR  REMAINING  SCENARIOS 


Scenario  4 (cont'd) 


73.5000 

0.0000 

0.0070 

0.0000 

74.0000 

0.0000 

0.0000 

0.0000 

74.5000 

0.0000 

-0.0600 

0.0000 

75.5000 

0.7000 

-0.0600 

0.0000 

76 . 5000 

0.0000 

-0.0600 

0.0000 

77.5000 

-0.7000 

-0.0600 

0.0000 

93.0000 

0.0000 

0.0000 

0.0000 

94.0000 

-0.5200 

0.0300 

0.0000 

95.0000 

0.0000 

0.0300 

0.0000 

96.0000 

0.5200 

0.0300 

0.0000 

110.0000 

0.0000 

0.0000 

0.0000 

111.0000 

0.5200 

-0.0300 

0.0000 

112.0000 

0.0000 

-0.0300 

0.0000 

113.0000 

-0.5200 

-0.0300 

0.0000 

120.0000 

0.0000 

0.0000 

0.0000 

Scenario  5,  AAl : 

42  discrete 

data  points 

2.0000 

0.0000 

0.0000 

0.0000 

3.0000 

0,4360 

0.0000 

0.0000 

4.0000 

0.0000 

0.0000 

0.0000 

5.0000 

-0.4360 

0.0000 

0.0000 

6.0000 

0.0000 

0.0000 

0.0000 

8.000C 

0.0000 

0.0610 

0,0000 

9.0000 

0,0000 

0.0000 

0.0000 

1 0.0000 

0.0000 

0.0000  - 

•10.0000 

11.0000 

0.1130 

0.0000  - 

•10.0000 

12.0000 

0.1130 

0.0000 

0.0000 

13.0000 

0.0000 

0.0000 

0.0000 

14.0000 

-0.1130 

0.0000 

0.0000 

15.0000 

-0.1130 

-0.0310 

0.0000 

16.0000 

-0.0700 

-0.0310 

10.0000 

17.0000 

0.0000 

0.0000 

10.0000 

21.0000 

0.0000 

0.0000 

0.0000 

23.0000 

0.0000 

0.0000 

2.0000 

29.0000 

0.0000 

0.0000 

0.0000 

30.0000 

0.0700 

0.0000 

0.0000 

31.0000 

0.0000 

0.0000 

0.0000 

32.0000 

0.3320 

0.0000 

0.0000 

33.0000 

0.0000 

0.0000 

0.0000 

34.0000 

0.0000 

0.0700 

0.0000 

35.0000 

-0.3320 

0.0000 

0.0000 

36.0000 

-0.0870 

0.0000 

0.0000 

37.0000 

0.0000 

-0,0700 

0.0000 

45.0000 

0.0000 

0.0000 

0.0000 

46.0000 

0.0870 

0.0000 

0.0000 

47,0000 

0.3140 

0.0000 

0.0000 

48.0000 

0,0000 

0.0000 

0,0000 

49.0000 

-0.3140 

0.0000 

0.0000 

i 
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DISCRETE  DATA  FOR  REMAINING  SCENARIOS 


Scenario  5 (cont'd) 


51 ♦ 0000 

0.0000 

0.0000 

0,0000 

52.0000 

0.0000 

0.0000 

-2.0000 

53.0000 

0.0000 

0.0700 

-2.0000 

56.0000 

0.0000 

0.0000 

0.0000 

57.0000 

0.3140 

0.0000 

0.0000 

58.0000 

-0.3140 

0.0000 

0,0000 

59.0000 

-0.3490 

0.0000 

0.0000 

60.0000 

0.1220 

-0.0350 

0.0000 

62.0000 

0.1220 

0.0000 

0.0000 

63.0000 

0.0000 

0.1400 

0.0000 

65.0000 

0.0000 

0.0000 

0.0000 

Scenario  6,  AA2: 

l6  discrete  data  points 

1.0000 

0.0000 

0.0000 

0.0000 

2.0000 

0.5410 

0.0000 

0,0000 

2.5000 

0.5410 

0.0350 

0.0000 

4.0000 

-0.5410 

0.0350 

0.0000 

6.0000 

-0.1400 

0.0350 

0.0000 

8.0000 

0.1400 

0.0350 

0.0000 

9.0000 

0.3500 

0.0350 

-24.3000 

10.0000 

0.0520 

0.0350 

-24,3000 

12.0000 

0.0000 

0.0350 

0.0000 

13.0000 

-0.4010 

0.0000 

0.0000 

14.0000 

-3.7520 

0.0000 

0.0000 

14.5000 

0.0000 

—3 • 6650 

0.0000 

15.0000 

0.0000 

3.6650 

0.0000 

16.0000 

3.7520 

0.0000 

24.3000 

17.0000 

1.0470 

0.0000 

24.3000 

18.0000 

0.0000 

0.0000 

0.0000 

Scenario  7 > M3: 

8 discrete 

data  points 

3.0000 

0.0000 

0.0000 

0,0000 

4.0000 

0.5240 

0.0210 

0.0000 

5.0000 

0.0000 

0.0210 

0.0000 

6.0000 

-0.5240 

0.0210 

0.0000 

7.0000 

-0.9600 

0.0210 

0.0000 

8.0000 

0.0000 

0.0210 

-25.0000 

9.0000 

0.9600 

0.3190 

-25.0000 

11.0000 

0.0000 

0.0000 

0.0000 
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APPENDIX  A 

NUMERICAL  INTEGRATION  TECHNIQUE 


In  general,  numerical  approximations  are  used  whenever  an  integration 
is  desired.  The  numerical  technique  used  in  this  software  is  given  below: 


z(t)  = z(t  - At)  + At/2  [3z(t)  - z(t  - At)] 


where 


z(t) 
z(t  - At) 
z(t) 
z(t  -At) 
At 


state  at  time  t sec 
state  at  time  t — At  sec 
dz/dt  at  time  t sec 
dz/dt  at  time  t — At  sec 
integration  interval  ( in  sec ) 


For  these  subroutines,  an  integration  interval  of  0.1  seconds  was  used 
(At  = 0.1). 

This  method  is  implemented  within  the  subroutines  by  saving  the  pre- 
vious values  of  z(t)  and  z(t)  in  the  arrays  z(t  - At)  and  z(t  - At).  In 
FORTRAN  variable  name  notation,  the  current  values  of  the  states  and  their 
derivatives  are  stored  in  the  Z and  ZD  arrays,  while  the  past  values 
(at  t — At)  of  the  states  and  their  derivatives  are  stored  in  the  ZO  and 
ZDO  arrays.  Thus,  the  FORTRAN  expression  which  implements  this  numerical 
approximation  is: 

Z « ZO  + DELT2*(3.0*ZD  - ZDO) 
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APPENDIX  B 


NOTES  OH  SOFTWARE  USE 


1 . To  execute  the  scenario  generation  package: 

a.  Place  appropriate  PROGRAM  card  in  front  of  NEWMAN  mainline 

routine.  ^ 

b.  The  following  routines  must  be  linked  to  execute  scenario 
operation: 


mainline 
SENGEN  ' 
DATGEN 
READAT 
HQDTU 
ERRMES 
OUTR  \ 
GPDPLT  ( 


(File  name  NEWMAN) 
(File  name  SENLIB) 
(File  name  PLTLIB) 


2.  To  execute  drive  logic  simulation/measure  generation  package: 


a.  Place  appropriate  PROGRAM  card  in  front  of  NEWMAN  mainline 
routine.  This  PROGRAM  card  should  contain  the  name  of  each 
tape  to  be  used/  since  each  set  of  scenario  data  resides 

on  a separate  tape  (conceptually,  at  least),  a tape  name 
should  appear  on  the  PROGRAM  card  for  each  set  of  scenario 
data  to  be  processed  through  the  drive  logic. 

b.  The  following  routines  must  be  linked  to  execute  drive 
logic  — sensory  apparatus  simulation 


ROLSWA 
or  PITSUR 
or  YAWHEV 
THRES  l 
THRESD  f 


(File  name  ROLSWA) 

(File  name  PITSUR) 

(File  name  YAWHEV) 

(File  name  DLLIB,  only  necessary 
if  ROLSWA  is  used) 


I 
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